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CAPÍTULO I 





1.1.- CONCEPTOS BÁSICOS 


1.1.1.- INTRODUCCIÓN 


Un paso previo al estudio de los métodos existentes para analizar la respuesta de un 
sistema estructural ante solicitaciones sísmicas, es la comprensión de las características de 
la respuesta de la estructura ante una solicitación dinámica de tipo general. 

El comportamiento dinámico difiere del estático en dos aspectos fundamentales. El 
primero de ellos se refiere a la dependencia con el tiempo de la solicitación, lo que lleva 
consigo la determinación de la respuesta en cada uno de los intervalos de tiempo de 
actuación de la carga. 

Sin embargo, la diferencia más significativa entre el comportamiento estático y 
dinámico de una estructura radica en la aparición, durante el movimiento, de fuerzas de 
inercia relacionadas con la masa del sistema. El juego entre la energía cinética de las masas 
y la energía potencial almacenada en la estructura es el factor que controla el movimiento y 
el cociente entre rigidez y masa está relacionado con una magnitud característica de los 
sistemas lineales: la frecuencia de oscilación. 

Un ejemplo clásico de la importancia que puede tener esta combinación de fuerzas elásticas 
y de inercia es el colapso producido en un puente de acero en Broughton (1831), por las 
vibraciones inducidas por un grupo de 60 soldados al atravesarlo marcando el paso; 
consecuencia de este accidente fue la colocación de carteles informativos (figura 1.1), 
indicando la conveniencia de romper el paso de las tropas cuando hubieran de atravesar un 
puente. 

En este ejemplo, los soldados aplican sobre el puente una fuerza con una determinada 
frecuencia que produce el movimiento del sistema; consecuencia de este movimiento y 
siguiendo la 2 a Ley de Newton, es la aparición de unas fuerzas de inercia proporcionales a 
la aceleración, estas fuerzas alejan al puente de su posición de equilibrio; por otra parte, 
entran en acción las fuerzas elásticas que deberán equilibrar no sólo el sistema de 
solicitaciones exteriores sino también a las fuerzas de inercia, devolviendo al puente a su 
posición de equilibrio. Resultando, en consecuencia, una oscilación del sistema alrededor 
de la citada posición con una frecuencia que depende de las características inerciales y 
elásticas del sistema; esta tendencia a la vibración se agrava, como se verá más adelante, 
cuando la frecuencia de oscilación del sistema coincide con la frecuencia de la excitación, 
como ocurrió en esta desafortunada situación. 
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Fig. 1.1 “Albert Bridge”Londres 

1 . 1 . 2 .- MOVIMIENTO PERIÓDICO 

OSCILACION O VIBRACION: Es el fenómeno que caracteriza el cambio cíclico 
de la posición de un sistema alrededor de una referencia a intervalos regulares de tiempo. 

Si se designa por x(t) al desplazamiento de un sistema, en el instante t; un 
movimiento vibratorio viene caracterizado mediante la ecuación: 

x(t)=x(t + T) (1.1-1) 

en la figura 1.2 se representa el desplazamiento, en función del tiempo, provocado en la 
cimentación soporte de una turbina de gas. 



Figura 1.2 
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PERIODO: Es el tiempo, T, que tarda una partícula del sistema en realizar una 
vibración completa, es habitual medirlo en segundos. 

En general si x(t) es una función periódica de período T, su desarrollo en serie de 
Fourier es . . . . . 
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Las funciones 


a n = eos ncot 
bn = sen ncot 


se denominan armónicos de orden n de la función periódica x(t), siendo su período T/n. 


1.1.3.- MOVIMIENTO ARMÓNICO 

La forma más simple de un movimiento periódico recibe el nombre de Movimiento 
Armónico, siendo la ecuación que describe el citado movimiento: 

x(t)=Asencot (Ll-2) 

donde: 

t representa el tiempo 

A valor máximo del desplazamiento, se designa por amplitud 
co frecuencia circular o pulsación, se expresa en radianes por segundo (rd/seg). 
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Se define la FRECUENCIA (f), como el número de oscilaciones completas realizadas 
en la unidad de tiempo, siendo la inversa del período, T; es decir: 



o> 

271 


(1.1-3) 


se mide en ciclos por segundo (HERTZ) 

El movimiento de la masa m provocado por el mecanismo indicado en la Figura 1.3, es 
un ejemplo de movimiento armónico. 

Este sistema, típico de las máquinas de coser, consta de una manivela de radio A, unida 
en uno de sus extremos a un punto O alrededor del cual gira; el otro extremo P desliza 
sobre una biela de guía vertical R. Al girar la manivela con una velocidad angular oo, el 
punto S y por tanto la masa m, se desplaza una cantidad que viene dada por: 


siendo la velocidad: 


x = AsenO = Asencot 
\ = (OA COSCDt 


y la aceleración: 

x=- © 2 Aseno)t = -oo 2 x 



Fig.1.3 
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1.1.4.- REPRESENTACION DE UN MOVIMIENTO VIBRATORIO ARMÓNICO 

El movimiento vibratorio armónico, suele representarse mediante la proyección, sobre 
una línea recta, de un vector OP, de módulo A que gira a una velocidad angular constante co 



Fig. 1.4 
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En la Figura 1.4, se representan las proyecciones horizontal y vertical del vector OP, al 
que se ha designado por X: 

• Proyección horizontal 

x(t) = A eos cot (1.1-4) 

® Proyección vertical 

x(t) = A sen ©t (1.1-5) 

Tomando la proyección horizontal, pueden presentarse las siguientes situaciones: 

(1): Si para t = 0, el punto P se encuentra en Q, el punto M proyección del P sobre el 
eje X, tendrá por abscisa: 

x = Acoscot 



(2): Si para t = 0, el punto P se encuentra en Q, entonces 

0 = ©t + 9 0 

siendo 0 O el valor de 0 para t = 0; y la proyección del punto P sobre el eje X, tomará 
el valor 

x = Acos (©t+0 o ) 
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(3): Si para t = O, el punto P se encuentra en Q, entonces 

x = Acos (cat — Bg ) 



El ángulo 0 O recibe el nombre de FASE INICIAL, o ángulo con el que se inicia el 
movimiento 

La velocidad y aceleración del movimiento armónico, se obtendrán diferenciando las 
ecuaciones (1.1-4) ó (1.1-5); así para un desplazamiento: 


x(t) = A eos cot 
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se obtiene 


x(t)=- A sencot = 


dx 

dt 


(1.1-6) 


x(t)=- Acó 2 cosco t = 



(1.1-7) 


A la vista de lo anterior, puede observarse que la velocidad y aceleración son también 
armónicas, con la misma frecuencia y amplitudes, Acó y Acó 2 , respectivamente; estando 
desfasadas, respecto al desplazamiento en n / 2 y n radianes. 


n 

x(t) = Acocos(cot H — ) 

2 

x(t) = Acó 2 eos (cot + 7t) 

Una representación de las relaciones (1.1-4), (1.1-6) y (1.1-7), se obtiene mediante el 
diagrama de ARGAND, indicado en la Figura 1.5. 



Fig. 1.5 
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Como puede observarse, la aceleración tiene signo opuesto al desplazamiento; si no 
fuese así, el movimiento vibratorio no se produciría; debe tenerse en cuenta que una 
aceleración negativa es una deceleración, esto es, el punto que oscila tiende a la posición de 
equilibrio. 

En la Figura 1.6, se indica la variación, con oot, del desplazamiento 


velocidad 
y aceleración 


u = u 0 sen cot 
ú = u 0 cocos oot 


ü = u 0 co 2 sen cot 



Fig.1.6 
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Un movimiento armónico admite también una representación mediante números 
complejos. En efecto, cualquier vector, X, en el plano XY puede expresarse mediante un 
número complejo 


siendo i, la unidad imaginaria 


X - a + ib 


i = 

y a, b las componentes real e imaginaria de x que representan sus proyecciones sobre los 
ejes x e y 



Si se designa por A al módulo de X y por cp a su argumento, puede escribirse: 

X = Ae 1 ’' = A(cos(p + isencp) (1-1-8) 

A = yja 2 + b 2 

b 

cp = arctg — 
a 


Diferenciando (1.1-8) se obtiene 


— = — (Ae i<Bt ) = iooAe ia> ‘ = iroX 
dt dt 

= — (ia>Ae i<Dt ) - -o 2 Ae '® 1 = -G) 2 X 
dt 2 dt 
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Si el desplazamiento viene representado como 

I m [Ae ,<II, ] - Asen ©t 

entonces 

velocidad = I ni [i©A e ,t0 ' ] = © A sen (©t + ^) 
aceleración = I m [-© 2 Ae tot ] = © 2 A sen (cot + ti) 

Es decir: “el movimiento es la proyección sobre el eje imaginario de un vector complejo 
de módulo A 

O bien, sí: 

desplazamiento = R c [Ae ,< “'] = A eos cot 

resulta que: 

velocidad = RjiwAe 1 ®'] = -©Asen ©t = ©Acos(©t +^) 
aceleración = R e [- © 2 Ae ,fflt ] = -© 2 A eos ©t - © 2 A cos(©t + 7t) 


1.1.5.- RELACION ENTRE LA FASE Y LAS CONDICIONES INICIALES 
Si la ecuación que caracteriza al movimiento puede describirse mediante la relación 

y(t) = A © sen( ©t + <j> ) 

y si se designa por yo e yo al desplazamiento y velocidad en el instante inicial, t = 0, se 
tiene que al ser 


y(t) = A © cos( ©t + <f> ) 


en t = 0, puede escribirse 


y o = A sen <|> 


y 0 = A © eos <j> 
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por lo que el ángulo de fase será 


y 0 ro 

<j) = arctag 

y« 


y la amplitud del desplazamiento, A, se obtendrá al tener en cuenta que 


y 0 


y « 12 


Aro 


= 1 


es decir 


A = ,y; + ^ 
V © 



o 


1.1.6.- MODELOS PARA SISTEMAS DINÁMICOS 

Se denominan coordenadas de una estructura, al número de parámetros que es preciso 
especificar para determinar su posición en cada instante de tiempo; y grados de libertad 
(g.d.l) a las coordenadas que resultan después de descontar las condiciones de enlace. Por 
las razones anteriormente expuestas, en el cálculo dinámico los g.d.l. a considerar serán 
aquellos directamente relacionados con la distribución de masa del sistema ; así, una 
estructura con una distribución continua de masa, tendrá un número infinito de g.d.l. 
Afortunadamente, en la mayoría de los casos, el movimiento del sistema estructural es tal 
que sólo es necesario especificar los N g.d.l. más significativos. A este tipo de sistemas se 
les denomina sistemas con N g.d.l.; utilizándose también, para el caso lineal, el apelativo 
oscilador lineal múltiple (o.l.m.). 
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La respuesta dinámica de algunos sistemas estructurales puede analizarse con la 
consideración de un único g.d.l.; es decir, sólo es necesario especificar una coordenada para 
describir su movimiento, lo cual no implica que no puedan producirse otros movimientos, 
sino que estos son despreciables frente al asignado a la precitada coordenada. 

En la figura 1.7. (a) se representa un modelo de una plataforma marina, en el que la 
pesada plataforma superior se ha reemplazado por una masa m, mientras que la estructura 
soporte proporciona la rigidez lateral suficiente contra el viento y oleaje; por lo que puede 
considerarse como un muelle de constante elástica k. 




y 


k 


-AA/V 

m 


FO) 


(a) 



Fig.1.7 
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De una forma análoga el puente representado en la Figura 1.7. (b), puede modelizarse 
mediante un sistema con g.d.l., que represente el movimiento vertical que sobre el centro 
del mismo provoca el paso de un vehículo. 

La Figura 1.8. (a) representa una construcción de una sola planta sometida a una 
explosión; un modelo de la estructura se indica en la Figura 1.8.(b), dicho modelo supone a 
los pilares flexibles a las deformaciones laterales y sin masa; concentrando la masa total del 
sistema en un dintel rígido. 



Fig.1.8 


En todos los casos citados, la solicitación exterior varía con el tiempo y los movimientos 
más representativos se reducen a uno sólo. 

En resumen, para el oscilador lineal, se supondrá la existencia de una sola masa cuya 
situación queda determinada por una variable de posición, y(t) , en cada instante t. Las 
derivadas respecto al tiempo de esta variable se indican mediante puntos; así, y(t) 
representa la velocidad e y(t) la aceleración 
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1.1.7.- ECUACIÓN DEL MOVIMIENTO 

El equilibrio de una partícula de masa m puede escribirse como: 

my = £F (1.1-9) 

donde: 

m : es la masa 

y : es la aceleración de la partícula 
F : las fuerzas que actúan sobre ella 

La expresión (1.1-9) puede interpretarse como el ESTADO GENERAL DE 
EQUILIBRIO DE UN SISTEMA DE RESULTANTE NULA, en cuyo caso el término, 
-m y , es la denominada FUERZA DE INERCIA. Siendo la presencia de estas fuerzas de 
inercia; esto es, la MASA, lo que marca la diferencia entre los problemas estáticos y 
dinámicos. 

Es usual considerar varios tipos de fuerzas F, y entre ellas: 

a. - Solicitaciones que actúan sobre la partícula f(t) 

b. - Reacciones elásticas -ky , donde k es una constante denominada 
COEFICIENTE DE RIGIDEZ 

c. - Reacciones disipativas del tipo -cy , siendo c el denominado 
COEFICIENTE DE AMORTIGUAMIENTO 

La Figura 1.9. (a) recoge el modelo típico de un sistema con 1 g.d.l 



( 


(a) 


Fig.1.9 
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Dicho modelo está formado por una masa, m, unida a un soporte fijo por dos elementos 
en paralelo. 

® Las características elásticas de la estructura se suponen representadas por la rigidez k 
de un muelle lineal de forma que para conseguir un desplazamiento relativo entre los 
extremos del mismo, es preciso actuar con una carga 

F k = ky(t) (U-10) 

• Se introduce en el modelo un amortiguador viscoso de constante c, tal que una 
velocidad relativa entre sus extremos, y(t) , sólo se consigue por aplicación de una 
fuerza 


F k =cy(t) (1.1-11) 

Este amortiguador es la imagen de la disipación de energía observada en los sistemas 
reales; aunque la forma de representación está más relacionada con la sencillez introducida 
en el modelo que con la realidad física. 

Otros tipos de amortiguamiento tienen un tratamiento matemático más complejo, que 
posteriormente será planteado en otros apartados. 

En la Figura 1.9.(b) se representan las fuerzas elásticas Fk, de disipación F c y de inercia 
Fi que actúan sobre el sistema excitado por una solicitación F(t) variable con el tiempo. 

La respuesta y(t) del sistema se obtendrá resolviendo la ecuación del movimiento cuya 
escritura se consigue al poner el equilibrio de fuerzas actuante sobre el sistema. Así: 

F r +F c +F k =F(t) 


donde 

F, = m y(t) 

Es importante hacer notar que el cumplimiento lineal de (1.1-10) y (1.1-1 1) hace lineal 
también la ecuación diferencial de equilibrio obtenida al sustituir ambas expresiones en 
(1.1-9) 


my( t )+ c y( t ) + k y( t ) =F ( t ) (li- 12 ) 

Es usual, reservar nombres especiales para casos particulares de la ecuación (1.1-12): 


VIBRACION LIBRE 
F(t) ~ 0 


SIN AMORTIGUAMIENTO : c = 0 
AMORTIGUADAS : c* 0 
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VIBRACION FORZADA SIN AMORTIGUAMIENTO : c = O 
F(t) * O AMORTIGUADAS : c *0 


1 . 2 .- VIBRACIONES LIBRES 

Se analiza en este apartado el movimiento vibratorio de sistemas sometidos a fuerzas 
inherentes a la propia estructura en ausencia de solicitaciones exteriores; es decir, cuando 
F(t) = 0 


1.2.1.- VIBRACIONES LIBRES SIN AMORTIGUAMIENTO 

En este caso, y al ser c = 0 , el sistema consta de una masa y un muelle tal como se 
indica en la Figura 1.10 



k§ 


1 


m 


T 


w 


k 


(5+yp 

— - 


m 



Fig.1.10 


La figura 1.10-(I), representa un muelle de masa despreciable, longitud L y 
constante k (razón de la fuerza al alargamiento). 
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Si se sitúa una masa, m, en el extremo A (Posición Inicial), Figura Í.IO-(II), se 
perturba el equilibrio del sistema, estableciéndose en una nueva posición A’ obtenida como 
consecuencia de un alargamiento 8 del muelle de forma que su longitud pasa a ser L + 8; 
esta nueva posición se habrá obtenido al considerar, por el Principio de Acción y Reacción, 
que el muelle ejercerá una fuerza hacia arriba igual y de signo contrario al peso del cuerpo, 
es decir, 

-kS-t-o) = 0 

o lo que es igual 


w = k 8 


( 1 . 2 - 1 ) 


donde 8 es la deformación del muelle en su posición de equilibrio estático y se denomina 
FLECHA ESTÁTICA. 

Si a partir de esta posición, y como consecuencia de algún desplazamiento, y(0), ó/y 
velocidad, y(0) , la masa m, se desplaza una cantidad, y, la partícula oscilará alrededor de 
la posición de equilibrio estático. 

Aplicando ahora la Segunda Ley de Newton: “ Cuando una fuerza actúa sobre un 
sistema produce en él un movimiento uniformemente acelerado, cuya aceleración es 
proporcional a la fuerza y de la misma dirección y sentido ”, puede escribirse 

Por otra parte, la resultante de las fuerzas aplicadas sobre la masa m será: 

2 F = w - k(8 + y) 

Igualando las dos últimas expresiones, y teniendo en cuenta (1.2-1) puede escribirse 

my + ky=0 (1.2-2) 

Obsérvese que al tomar como referencia para, y, la posición de equilibrio estático, se 
elimina de la ecuación del movimiento: el Peso, w, y la Fuerza estática sobre el muelle, kS 

Dividiendo por m los dos miembros de (1.2-2), se obtiene: 

y + o 2 y = 0 (1.2-3) 


donde co 2 = K/m, y o es la denominada FRECUENCIA NATURAL DEL SISTEMA, que 
representa una medida entre los parámetros ELÁSTICO, k, e INERCIAL, m, del sistema. 
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La solución de (1.2-3) podrá encontrarse con soluciones particulares del tipo: 

y =Ce ri 

y = C r e rt (1.2-4) 

y =Cr 2 e rt 


Sustituyendo (1.2-4) en (1.2-3) se obtiene la ecuación característica 


cuya solución es: 
por lo que 


o bien 


r 2 + cd 2 = 0 


= ± o) V-T = ± i 


ico 


_ ícot 


y = C, e'" 0 ' + C 2 e“ 
y = Asencot+ B coscot 

siendo A y B constantes que dependerán de las condiciones iniciales; si éstas son 


se obtiene para (1.2-3) 


y(t = O) = y 0 

y(t = o) = y 0 


y(t) = y 0 eos cot + 


V (0y 


sen cot 


Esto es, el sistema vibrará tal y como se indica en la Figura 1.11. 

El sistema inicialmente en su posición de equilibrio (1) se mueve hacia la derecha hasta 
alcanzar el máximo desplazamiento, posición (2), a partir de la cual comienza a disminuir 
hasta alcanzar de nuevo la posición de equilibrio (3), punto de comienzo del movimiento 
del sistema hacia la izquierda hasta alcanzar la posición (4) en la que los desplazamientos 
comienzan a disminuir hasta alcanzar la posición de equilibrio (5), repitiéndose de nuevo el 
ciclo completo. 
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Fig. 1.11 Deformada del sistema en las posiciones 1, 2, 3, 4 y 5 


Teniendo en cuenta lo anterior, es posible establecer un BALANCE DE ENERGÍA, lo 
que se consigue al multiplicar los dos términos de la ecuación (1.2-2) por la velocidad, y . 


o bien: 


my y + kyy = 0 



+ k y + 


dy 

dt 


= 0 


e integrando respecto del tiempo 



^dt + ík y + — 
dt 3 dt 


dt 


= 0 
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se obtiene la expresión del Principio de Conservación de la Energía, es decir: 

1 .2 1 . 2 ^ 

— rn y + — k y = E = cte. 

2 J 2 J 

donde, la Energía almacenada en el muelle producida por la deformación elástica viene 
representada por el término 

y la almacenada en la masa en virtud de la velocidad 

E = — m y 2 
c 2 J 

Como para sistemas conservativos la energía ha de permanecer constante: 

• Cuando la velocidad del sistema es nula, la energía cinética, E c , es nula, alcanzando 
la energía potencial su valor máximo, E pmx 

Ep mx = E 

• De una forma análoga a desplazamientos nulos la energía potencial desaparece 
tomando la Energía Cinética su valor máximo, E pmx 


Ecmx = E 


Es decir, si la energía total del sistema permanece constante, entonces la Energía cinética 
máxima es igual a la potencial máxima. 

E( mx Ep mx 

Con lo cual se pone de manifiesto que el movimiento del sistema es un continuo 
intercambio de energía cinética y potencial entre la masa y el muelle. 
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1.2.2.- PRESENCIA DE AMORTIGUAMIENTO. DISIPACIÓN DE ENERGÍA 

El proceso por el cual, la vibración disminuye en amplitud, recibe el nombre de 
amortiguamiento. En este proceso la energía del sistema se disipa como fricción o calor. 
Por lo general un sistema presenta siempre una cierta amortiguación en su movimiento que 
a veces es despreciable; sin embargo, en general, usualmente es inevitable su 
consideración. Entre los distintos tipos de amortiguamiento, cabe señalar: 

L- AMORTIGUACION FLUIDA 

Puede ser viscosa o turbulenta: 

• VISCOSA: la fuerza de amortiguación es proporcional a la velocidad 

• TURBULENTA: la fuerza es proporcional al cuadrado de la velocidad 

2. - AMORTIGUACION DE FRICCIÓN (COULOMB) 

La fuerza de amortiguamiento es constante 

3. - AMORTIGUACION HISTERETICA 

También denominada amortiguación estructural, está producida por la fricción 

interna o histéresis cuando se deforma el sistema. 

1.2.2.1- Amortiguamiento viscoso 

Un medio fluido presenta una resistencia al movimiento de un sistema, produciéndose 
una disipación de energía que depende de factores tales como la forma, dimensiones y 
frecuencia de vibración del fluido. 

Un mecanismo de amortiguación usual es el amortiguamiento viscoso, en el cual la 
fuerza es proporcional a la velocidad; realmente esto sólo es válido para una amortiguación 
como la producida por el flujo laminar de un fluido viscoso. Sin embargo, otras formas de 
amortiguamiento se aproximan a este último tipo si las fuerzas de disipación son pequeñas. 

La utilización del amortiguamiento viscoso, tiene la ventaja de linealizar la ecuación del 
movimiento. 

El amortiguador viscoso, se representa esquemáticamente por un pistón y la constante c, 
tal y como se indica en la Figura 1.12. 

c 



Fig. 1.12 
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Para el caso de vibraciones libres amortiguadas (amortiguamiento viscoso), el equilibrio 
del sistema (Figura 1.13) implica: 


my + cy + ky = 0 (1-2-5) 

que suele escribirse en la forma: 

y + 2C<oy + (Q 2 y = 0 (1.2-6) 



Fig.1.13 


Donde se ha hecho 

2 k 

<» =— (1.2-7) 

m 

r = - c 

V4km 

siendo: C, el ÍNDICE DE AMORTIGUAMIENTO 

La ecuación (1.2-6) es una ecuación diferencial ordinaria cuya solución es del tipo: 

y(t) = Ce rt (1.2-8) 

sustituyendo (1.2-8) en (1.2-6), se obtiene: 

r 2 + 2^cor + co 2 = 0 (1-2-9) 

es decir 


r - ® £~ C 
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Debido a linealidad del problema la solución general será: 

y(t) = e _í “* ^Ae'^* <ot + Be"^®' 


(1.2-10) 


dependiendo las constantes A y B de las condiciones que desencadenan el movimiento. 
Obsérvese que el primer factor de la ecuación, é^ at , es modulador del segundo y disminuye 
con t. 


La ecuación (1.2-10) adopta formas distintas según sean reales o imaginarias las raíces n 
y r 2 de la ecuación (1.2-0); lo que se corresponde físicamente con una amortiguación 
grande (¿¡ > 1), ó con una amortiguación débil (C, < 1); existiendo entre estos dos casos 
extremos la amortiguación crítica correspondiente a £ = 1. 

• Sistema Sobreamortiguado (£ > 1) 

En este caso, las dos soluciones ri y r 2 de la ecuación (1.2-9), son reales 

r 2 =[-^-^~í\ o) 


por lo que la ecuación (1.2-10) puede escribirse en la forma 


, . , _ i hot I -C-VC- 1 1 ®>t 

y(t) — C, e L J + C 2 e L 

determinándose Ci y C 2 a partir de las condiciones iniciales; si éstas son: 

y(í = 0 ) = y 0 


y(t = 0 ) = 


c,= 

C 2 = 


y 0 r 2 -y 0 

y o fí-y 0 
r !~ r 2 


se obtiene: 
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o bien 



Es decir, el sistema al iniciar el movimiento se verá amortiguado aproximándose 
exponencialmente a la posición de equilibrio, que teóricamente nunca llegará a alcanzar, no 
produciéndose, en este caso, oscilación alrededor de la citada posición. 



• Sistema Críticamente amortiguado (£ = 1) 

Si C, = 1, las dos raíces de la ecuación (1.2-9) son iguales a, -oo, la ecuación del 
movimiento será: 

y(t) = e“'[C 1 t + Cj 

La aplicación de las condiciones iniciales 

y(t = 0) = y 0 
y(t = 0) = y 0 

proporciona, en este caso, para Ci y C 2 

c, =y 0 +®y 0 
c 2 =y 0 
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siendo la ecuación del movimiento 


y(t) = e “'[y 0 +(y 0 +®y 0 )t] 

por lo que el sistema vuelve a la posición de equilibrio sin oscilar 



c = V4 k m 

este valor de la constante de amortiguamiento, recibe el nombre de 
AMORTIGUAMIENTO CRITICO y depende de las características inerciales y de rigidez 
del sistema 

C crit = V4km (1.2-11) 

Teniendo en cuenta la ecuación (1.2-7), puede escribirse: 

S = — (1.2-12) 

C crit 

por tanto, el índice de Amortiguamiento £ es una medida adimensional del coeficiente c 
de amortiguamiento del sistema. 

• Sistema Sub-amortiguado 

Si ^ < l,o lo que es igual 

« < C cri< 

se producen para r soluciones imaginarias conjugadas 


r = ©[-¿;± iyjl-C , 2 1 
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y la solución para y(t) será: 


y(t) [Ae it0l,t +Be _i “ d '] (1.2-13) 

donde co d es la frecuencia del sistema amortiguado 

c C (1.2-14) 

la ecuación (1.2-13) se puede escribir también, reordenando términos en la forma: 

y(t) = e” c<Dt [C, cosco,, t + C 2 sen to d t] (1.2-15) 

expresión que representa un movimiento vibratorio alrededor de la posición de equilibrio 
con una amplitud que disminuye con el tiempo y cuya frecuencia circular es coa , expresión 
(1.2-14). 


En este caso, el desplazamiento en el instante inicial, y(0), se ve modulado por una 
exponencial decreciente. Mientras que la cantidad entre corchetes se repite a intervalos: 


T 


D 


2n 


(O 


d 



Fig. 1.14(c) 
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De la ecuación (1.2-14) se deduce que la presencia del amortiguamiento en el sistema, 
disminuye su frecuencia natural de vibración, aumentando por tanto su período. Una 
representación de la citada ecuación, es la indicada en la Figura 1.15, donde puede 
observarse que para valores pequeños de amortiguamiento, la diferencia entre coa y © es 
despreciable a efectos prácticos. 



Fig. 1.15 


En la Tabla 1.2.1 se indican algunos valores característicos del índice de 
Amortiguamiento: 


1 

C(%) 

Automóviles (absorción en el choque) 

10-50 

Suelo 

5 

Goma, Caucho 

4 

Vigas de hormigón 

4 

Vigas de hormigón pretensado 

1,5 

Estructuras de acero remachadas 

1,5 

Estructuras de acero soldadas 

0,5 

Paredes de ladrillo 

0,8 


Tabla 1.2.1 
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1.2.2.2- Decremento logarítmico 

Una cantidad de interés, indicativa del grado de modulación introducido por el 
amortiguamiento en los desplazamientos del sistema se obtiene al comparar dos 
elongaciones y(t) e y(t+To), separadas por el período T D . 


Según (1.2-15): 


y(0 = e OT D 

y(t + T„) 


(1.2-16) 


y 


8 = | n-^~rr = C» T D (1.2-17) 

y(t + T„) 

la expresión (1.2-17) define el denominado DECREMENTO LOGARÍTMICO §; pudiendo 
escribirse, teniendo en cuenta (1.2-15) y (1.2-14), en la forma 


8 = 2 7C 


(1.2-18) 


siendo, por tanto, 8 una medida del “ grado de decrecimiento de la oscilación ”. 

El mismo razonamiento puede utilizarse para j períodos; en efecto, si se designa por y¡ , 
y¡+j las amplitudes en los instantes t¡ y t¡ + j T D , respectivamente se obtendría: 


y¡ y(t,) 
y¡ +j y(t,+jT D ) 


(1.2-19) 


tomando logaritmos en (1.2-19) y teniendo en cuenta (1.2-17) puede escribirse: 


ln-^- = C(ojT D =j8 (1.2-20) 

y i+j 

La expresión anterior, junto con (1.2-18), proporciona: 

8 = - In -^- = 2 ti - t J ?- . (1.2-21) 

J y ¡+ j vi -c 2 

resultado que puede utilizarse para evaluar experimentalmente el amortiguamiento de un 
sistema, en la forma siguiente: 
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1- El sistema inicialmente en su posición de equilibrio es perturbado a una 
posición yo- 

2- Se registra la vibración producida en el sistema. 



3. - Se mide el tiempo utilizado para completar un ciclo; es decir, se determina 

el período Td- 

4. - Se mide y¡ , y¡+j . 

5. - Se calcula 



j y¡ + j 


6 - Se calcula C, 

Aunque el procedimiento es aplicable con j = 1, puede proporcionar resultados 
incorrectos para sistemas con un amortiguamiento pequeño. 


Por último, cabe indicar que como en las estructuras habituales el valor de C, es 
relativamente pequeño, está en general justificado simplificar la expresión del decremento 
logarítmico, tomando como definición para el mismo: 8 = 27t£ 



c ti. . 




Fig. 1.16 
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1.2.2.3- Consideraciones energéticas 


Con objeto de analizar los efectos del amortiguamiento desde otro punto de vista, es 
interesante reconsiderar la ecuación 


my + cy + ky = 0 

multiplicando por ye integrando respecto al tiempo 

T T T 

Jmy ydt+ Je ky dt = -Jcy 2 dt (1.2-22) 

0 0 o 


se obtiene la Energía Disipada por ciclo, Ed, donde 

T 

E d = Jcy 2 dt (1.2-23) 

0 


como c e y 2 son esencialmente cantidades positivas, el signo - indica que Ed es la Energía 
disipada; existiendo, por tanto, una continua pérdida de energía durante el movimiento. 

En el caso de amortiguamiento viscoso y suponiendo un movimiento armónico del tipo* 

y(t) = A sen at (1.2-24) 

la Energía Disipada, Ed, en un ciclo completo será: 



= n ca Á 1 


(1.2-25) 


y sí el sistema vibra a su frecuencia natural, a = co entonces 

E d = 7t(2 7t^©m)G)A 2 = 2nC,kA 2 (1.2.25a) 


En el apartado (1.3.2), se verá que un movimiento de estas características, corresponde a la respuesta de un 
sistema amortiguado sometido a una fuerza de excitación armónica de frecuencia a 
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Es decir, la Energía Disipada es proporcional al cuadrado de la amplitud, A, del 
movimiento; conviene observar que la citada energía depende de la frecuencia, no siendo 
constante para valores ñjos de la amplitud y del amortiguamiento. 

Por otra parte, al ser la Fuerza de disipación: 

F d = c y = c A a eos at (1 .2-26) 

y teniendo en cuenta (1.2-24); la relación entre Fd e y, es una elipse de ecuación 

(4-) 2 +(^V) 2 = ! (1.2-27) 

A ca A 



cuya área representa la Energía Disipada por ciclo. 

La relación entre la energía disipada en cada ciclo del movimiento y la energía total del 
sistema se la define como CAPACIDAD ESPECÍFICA DE AMORTIGUAMIENTO, \|/ 


V = 



Por otra parte, como la energía total del sistema, E, puede expresarse mediante la energía 
potencial máxima, o bien como la energía cinética máxima, resultará. 


E = E 


cmx 
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obteniéndose para vj/ 


y como 


siendo 



271 c 
m a 


— = 2 C© 

rn 



rn 


por lo que, para a igual a la frecuencia natural co del sistema 



Finalmente, teniendo en cuenta la definición de Decremento logarítmico, la capacidad de 
amortiguamiento para índices de amortiguamiento bajos será 

E 

vj/ = — « 4 ti ^ »í 2 8 (1.2-28) 

E 

lo que sugiere un procedimiento para calcular C,. 

La energía disipada por ciclo proporcionado por la ecuación (1.2-25), también es válida 
cuando se considera al amortiguador viscoso dispuesto en paralelo con un muelle de 
constante k; en efecto, la fuerza necesaria para provocar el movimiento será. 


F = ky + cy 

(1.2-29) 

y para un movimiento armónico de amplitud A y frecuencia a 


y(t) = Asenat 

(1.2-30) 

se obtendrá al sustituir en (1.2-29) 


F- ky = cAa eos at 

(1.2-31) 

o bien 


F-ky 

= eos at 

c Aa 

(1.2-32) 

y teniendo en cuenta 


y 

— = sen a t 

A 

(1.2-33) 



1-34 


la suma de los cuadrados de (1.2-32) y (1.2-33) permite escribir 


F-ky 

cAa 



= 1 


(1.2-34) 


• • r 2 

que es la ecuación de la elipse representada en la Figura 1.18 y que encierra un área rccA 
igual a la Energía Disipada, Ed. 



Fig. 1.18 

La Energía disipada en el ciclo completo será 

T T 

E d = |Fydt= J(ky + cy)ydt 
0 0 

y como ky es perpendicular a no producirá trabajo, por lo que 

T 

E d = Jcy 2 dt = Tiac A 2 

o 

Así pues, medir el área encerrada por la elipse es equivalente a calcular la energía 
disipada. 

La tangente a la elipse en el punto P(F = ky; y = A) es vertical, lo que sugiere una 
construcción gráfica sencilla para medir experimentalmente el valor de k y obtener para 
ello obsérvese que (1.2-34) también puede escribirse en la forma 

F = ky ica^A 2 -y 2 


(1.2-35) 
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y para F = O, se tiene que y = d 

-kd = ±ca -/A^-d 2 " 
k 2 d 2 =c 2 a (A 2 -d 2 ) 

(k 2 +c 2 +a 2 )d 2 = c 2 a 2 A 2 

es decir 

d _ ca 

A / c a 

k l/ ,+( T> 

Por otra parte, teniendo en cuenta que: 

cco_ eco _ c 
k m co m co a> 

la expresión (1.2-36) puede escribirse, para a igual a la frecuencia natural: 

d 2 ^ 

A V 1 + 4 C 2 

y para valores de ^ pequeños 

d - 2^ A 


(1.2-36) 


(1.2-37) 


(1.2-38) 


(1.2-39) 


expresión que sugiere de nuevo la posibilidad de medir £ por inspección de la curva de 
histéresis. 

En lo que se refiere a la medida de la constante k, esta puede obtenerse teniendo en 
cuenta que derivando (1.2-34), respecto a y 


F-ky 


(c a) 


dF 

dy 


+ y = o 


el máximo se obtendrá imponiendo la condición 


es decir 



k(F-ky) = c 2 a 2 ,+ y 
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por lo que 


kF = y(c 2 a 2 + k 2 ) 


y la pendiente secante, ms 



Fig. 1.19 


Recordando que para la frecuencia natural, (1.2-37) proporcionaba: 


C (O 

k 


2 ? 


m s puede escribirse 


m s =k(l + 4C 2 ) (1.2-41) 

En resumen, conociendo m s y d pueden obtenerse C, y k mediante las expresiones (1.2- 
39) y (1.2-41). 
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1.2.2.4- Amortiguamiento histerético 

La observación de la ecuación (1.2-25), permite concluir que la energía disipada por 
ciclo es independiente de la frecuencia o disminuye muy suavemente con ella. En la Figura 
1.20, se indica la variación de la Energía disipada con la frecuencia. 



Fig.1.20 

La experimentación indica una tendencia contraria a la aceptada en el planteamiento 
realizado para el amortiguamiento viscoso, razón por la que se acostumbra a suponer que 
“la energía perdida es independiente de la frecuencia”, lo que implica que la Fuerza de 
disipación debida al amortiguamiento es proporcional a la velocidad e inversamente 
proporcional a la frecuencia. 

El modelo de amortiguamiento en este caso es designado como “ interno , estructural o 
de histéresis”, ya que puede relacionarse con el ciclo de histéresis que presenta el diagrama 
tensión-deformación de un sistema con amortiguamiento que oscila alrededor de su 
posición de equilibrio estático, Figura 1.21, y para el que el área encerrada por el bucle de 
histéresis representa la energía disipada en un ciclo de vibración. 



Fig.1.21 
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Para ello, el producto ca se sustituye por una constante h, denominada CONSTANTE 
DE AMORTIGUAMIENTO HISTERETICO: 

h = c a = cte. (1.2-42) 

y como para la frecuencia natural del sistema 

c = 2^ © m 

la ecuación (1.2-42) puede escribirse 

h = 2 C, co 2 m = 2<¡¡— m 

m 


por tanto 


h = 2^k 


(1.2-43) 


y la energía Disipada, Ed, ecuación (1.2-25) 

E d = 7r h A 2 (1.2-44) 

es función únicamente de la constante, h y de la amplitud del movimiento, A; siendo 
directamente proporcional al cuadrado de esta última, lo que implica claramente la 
consideración de un problema no lineal. 

Recordando que para valores de C, pequeños, (1.2-39), proporcionaba 

d = 2£ A 

y teniendo en cuenta (1.2-43), se obtiene 


d__ h_ 

A ~ k 


(1.2-45) 


expresión que permite la identificación de h. 

En la hipótesis de movimiento armónico, el grado de modulación por ciclo introducido, 
por la consideración de amortiguamiento histerético, puede obtenerse mediante el balance 
de energía entre dos puntos, P y R, separados medio ciclo. 
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y el índice de Amortiguamiento equivalente a un amortiguamiento viscoso C, se obtendrá 
teniendo en cuenta (1.2.21) 


s * r nh 
8 h =27r¿> — 

es decir 


Normalmente no es posible separar la rigidez de una estructura de su amortiguamiento 
histerético, motivo por el cual en el modelo matemático que se utiliza para caracterizar el 
comportamiento del sistema, las dos características: rigidez y amortiguamiento reciben un 
tratamiento conjunto mediante lo que se denomina RIGIDEZ COMPLEJA. 

En el caso más general un movimiento armónico es del tipo 

y = A e iw ‘ (1.2-47) 


lo que permite trabajar simultáneamente con la proyección horizontal y vertical del mismo. 
Teniendo en cuenta (1 .2-47), la fuerza F vendrá dada por la expresión: 



F = ky + cy 

F = kAe im ' + iA(De imt 
F = (k + ia>c)Ae im ‘ 

(1.2-48) 

o lo que es igual 

F = (k + i co c) y 

(1.2-49) 

o bien 

F = k * y 


donde 

k* = k + koc 

(1.2-50) 

y como h = co c, resulta 

L 

k* = k + i h = k 1 + i — 
k 



donde la parte real de k* representa la rigidez y la parte imaginaria el amortiguamiento. 
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I.2.2.5.- Amortiguamiento de Coulomb 


Otro modelo de disipación de energía, corresponde a la amortiguación por fricción o 
amortiguamiento de Coulomb. 

En este caso la fuerza de amortiguamiento es constante en magnitud, y por tanto 
independiente de la amplitud y frecuencia, pero de sentido contrario a la dirección de la 
velocidad, por lo que el sentido de la fuerza cambiará cuando cambie la dirección del 
movimiento; motivo por el que es necesario plantear dos ecuaciones que caractericen el 
movimiento del sistema en cada una de las dos situaciones planteadas en la Figura 1.23. 



Siendo 

F d = -j^rNtagtf» (1.2-51) 

m 


donde (j> es el ángulo de fricción y N la fuerza normal. Valores típicos de <j> , se indican en 
la Tabla 1.2.2 


Material 

♦ 1 

Madera sobre madera 

20° 

Acero sobre acero 

10° 

Arena y Grava 

o 

O 

en 

Arcilla 

20° 
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Para el primer caso, Figura 1.23(b), la 2 a Ley de Newton permite escribir 

my + ky-Ntag<j> = 0 ; si y<0 (1.2-52) 

cuya solución es 


y = A, eos co t + B, sen o> t + 


N tag <(> 


(1.2-53) 


y para el segundo, Figura 1.23. (a) 

my + ky + Ntag<|> = 0 ; si y>0 

cuya solución es 


y = A 2 eos co t + B 2 sen co t 


N tag<|> 


(1.2-54) 


(1.2-55) 


donde las constantes Ai, Bi, A 2 y B 2 dependen de las condiciones iniciales en cada 
semiciclo sucesivo. 

Para el caso en el que las condiciones iniciales sean 

y(t = o) = y 0 

y(t = 0) = 0 


se tendría 


y« = a, + 

0 = B, co 


N tag <j) 


o lo que es igual 


A i = Yo ' 


B, = 0 


N tag <(► 


por lo que la ecuación del movimiento toma la forma 


, Ntagdr Ntagd> 

y = (y„ — 7 a1 ) eos» t+ — - — 


(1.2-56) 


solución válida para 0 < cot < tc ; cuando cot = ti 


, 7t . . 2Ntag<l). 

y(t = - ) = -(y« r-^) 

(o k 


(1.2-57) 
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71 


y(t=— ) = o 


00 


es decir, la amplitud inicial, yo , se ha reducido en una cantidad igual al doble del 
desplazamiento característico ^ ^ 


Las condiciones (1.2-57), son las iniciales para la siguiente mitad del ciclo, por lo que 
00 k k 


71 


y( t =-) = 0 = B 2 <D 


00 


lo que implica 


a 2 =y 0 - 


3 N tag (|> 


B 2 =0 


siendo entonces la ecuación del movimiento, para ti < cot < 2 n 

, 3 N tag ó. A Ntagó 

y = (y 0 2_l) eos g) t - - 


y al final del semiciclo, cot = 2 7t 


oo k 


y(t=— ) = o 

00 


(1.2-58) 


que representan las condiciones iniciales para el siguiente semiciclo. 


En la Figura 1.24 puede observarse que, en cada semiciclo, se reduce la amplitud en un 
2 N tag <b 

factor por lo que para un número de semiciclos, r , la amplitud del movimiento 


sera 


Yo 


2 N tag <j) 


cesando el movimiento del sistema al final del semiciclo para el cual la amplitud es menor 
que el desplazamiento característico, es decir 
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y 0 - r 


2 N tag (¡> 
k 


< Ntag(|) 
k 



Fig.1.24 


o lo que es igual 


y 0 


2 N tag (j> 
k 

2 N tag <¡) 
k 


Como en cada ciclo se reduce la amplitud en un valor igual a 

4 N tag <|> 
k 

el decremento será lineal 


y (t + T) = y (t ) 


4 N tag (j) 


(1.2-59) 


(1.2-60) 
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siendo entonces la pendiente de la recta moduladora 

-4 N tag <j> 


271 

© 


2N© tag(() 
ti k 


(1.2-61) 


El resultado anterior puede utilizarse para definir un amortiguamiento viscoso 
equivalente. De hecho la ecuación (1.2-16) puede escribirse 


y(t + T p ) 

y(t) 


- n -0»T D 


= 1 -£¡0)T d + 


(W ) 2 

2 ! 


+ ••• 


y como T d , para índices de amortiguamiento bajo es igual a 27t/co resulta 

y(t + T D ) = y(t)-27iGy(t) + ... (1.2-62) 


comparando (1.2-60) y (1.2-62), y haciendo y(t) = A; se obtiene 

2 ,^- 4 ^ 

K 

es decir 


4 N tag <j) 



donde A es la amplitud del movimiento 


(1.2-63) 


Por otra parte, según (1.2-25a), la energía disipada es 

E d = 2 7r A 2 k ^ eq = 4 A (N tag <j> ) (1.2-64) 

lo que puede interpretarse como el cuádruplo del trabajo realizado por la fuerza, N tag <¡>, en 
un cuarto de período. 
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1.3.- VIBRACIONES FORZADAS 

Normalmente, la vibración de un sistema se produce como consecuencia de la aplicación 
sobre el mismo de fuerzas que varían con el tiempo. La respuesta del sistema bajo este tipo 
de excitación se denomina vibración forzada, y corresponde a la solución de la ecuación 

m y(t) + c y (t) + k y (t) = F(t) (1.3-1) 

Como es habitual, la solución a la ecuación diferencial anterior, se obtendrá en la forma 

y(t) = y„(t) + y p (t) (1.3-2) 

donde, y 0 (t), es la solución de la ecuación diferencial homogénea, véase (1.2-15), 

y 0 (t) = e~ ?mt [C, eos <o d t + C 2 sen o d t] 

e, Yp 0)? I a solución particular que dependerá de la expresión de F(t). 


1.3.1. CARGA CONSTANTE 

El caso más sencillo se produce cuando 

F(t) = P = cte 



Fig.1.25 



1-47 


Por lo que en este caso 


y P (t) = 


p 

k 


y por tanto 


y 0 (t) = -r + e ?ra, [c,coso) d t + C 2 sen(D d t] (1.3-3) 

K. 


si se consideran como condiciones iniciales 


y(t) = o 
y(t)=o 


entonces 


c,=-£ 

k 


C, =- 


PC® 

k CD„ 


tomando, en este caso, la ecuación (1.3-3) la forma 


y(t) 


1 -e~^ mt (cos ® d t + p sen o d t) 

Vi-í 2 


(1.3-4) 


y para el caso en el que C = 0, se obtiene 


y(t) = — [l -coscot] 
k 


(1.3-5) 


Al cociente P/k, se le designa por el término FLECHA ESTÁTICA, y representa el 
desplazamiento que se produciría en el sistema al aplicar la carga P estáticamente 


(1.3-6) 


La representación gráfica de (1.3-5) se indica en la figura (1.26), y en ella puede 
observarse que la amplitud del desplazamiento, y mx , es 


y mx =2- = 25 st 
k 
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1.3.2.- RESPUESTA A CARGAS ARMONICAS 

Debido a que es posible expresar cualquier función periódica como suma de funciones 
armónicas; es importante considerar la respuesta de un sistema sometido a una solicitación 
de estas características. 


En este caso 

F(t) = P 0 eos <D 0 t 

y la solución particular, y p (t) será 

y p = A eos <¡o 0 1 + B sen <o 0 1 


(1.3-7) 


(1.3-8) 


que sustituida en (1.3-1), y de acuerdo con (1.3-7) proporciona para A y B los siguientes 
valores 


G ) 2 


(0, 


(1 “ i) + 4 ^ 


co 


Pq 1-P 2 

k (l-p 2 ) 2 +4^ 2 p 2 


2C^ 

Pq ^G) Pq 2CP 

k (1 _^ +4 ^ C0 1 k (1-P 2 ) 2 +4C 2 P 2 

GT G ) 2 


(1.3-9) 
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Por otra parte, si se hace 

g 

A = Y 0 eos <|> 0 tag <¡ > 0 = — 

A 

=> 

B = Y 0 sen <|> 0 Y 0 = Va 2 +B 2 

la ecuación (1.3-8) puede escribirse en la forma 


Donde 


y P (*) = Y, eos (oo 0 t -<►„) 


tag<J) 0 = 


2? 


1- 


®0 


CD 



(1.3-10) 

(1.3-11) 


(1.3-12) 


Ye- 


co 


( l _^|) 2 + 4 ; 2t0 ' 


co 


(1.3-13) 


siendo, 5 st , la FLECHA ESTATICA, igual a P/k. 

En resumen, la respuesta y(t) del sistema, se compone de dos términos que corresponden 
a: 


la solución homogénea que, como se vio en el apartado (1. 2.2.1), tiende a cero 
si t -» oo, por lo que yh(t) representa la componente transitoria del movimiento. 

La solución particular, con una amplitud Y 0 y un ángulo de fase <j>o, es la 
componente permanente del movimiento-, esta solución no está modulada por 
ninguna función exponencial, de manera que el movimiento perdura en el tiempo. 


Al cociente -2- se le denomina FACTOR DE AMPLIFICACION DINÁMICO 

8„, 



® 0\2 - "2 


a— fr+4c 


© 


2 

© 2 


(1.3-14) 
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Las figuras 1.27 y 1.28 son la representación gráfica de las expresiones (1.3-12) y (1.3- 
14); y en ellas puede observarse lo siguiente: 

• Tanto el ángulo de fase, (j> 0 , como el factor de amplificación D, dependen únicamente 
de la relación de frecuencias y del índice de amortiguamiento t,. 

• El amortiguamiento disminuye el factor de amplificación, D, para cualquier relación 
de frecuencias. En efecto, si se designa la relación de frecuencias por 


se obtendría 


P = 


®0 

G) 


d = [( i - p 2 > 2 + ( 2 ¿; P ) 2 ] 1/2 


por lo que D mx se produciría para una relación de frecuencias 

CD 


y sería igual a 


D„ 


1 


(1.3-15) 


• En el caso en el que p = 0, el factor de amplificación es igual a la unidad, y en 
consecuencia 



es decir, si la frecuencia de excitación, gi 0> es muy pequeña comparada con la frecuencia 
natural del sistema, ro, resulta que la amplitud del movimiento Y 0 , es igual a la flecha 
estática, o lo que es igual, el fenómeno dinámico es despreciable, estando controlado el 
movimiento del sistema por la rigidez del mismo. 


• Cuando p = 1, la frecuencia de excitación, co 0 coincide con la frecuencia natural del 
sistema, co, siendo en este caso 


D = 


2C 


en otras palabras, el factor de amplificación está controlado únicamente por el 
índice de amortiguamiento, independientemente de la masa o de la rigidez. 
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Así por ejemplo, para el caso en el que £, = 2%, se obtendría un factor de 
amplificación igual a 25, lo que supone una amplitud de desplazamiento 25 veces 
superior a la que se obtendría en el caso de aplicación estática de la carga. 

Esta amplificación máxima de la respuesta recibe el nombre de RESONANCIA, 
fenómeno que puede producir el colapso del sistema. 


La relación de amplitudes en la resonancia se denomina FACTOR DE CALIDAD, Q, 
término utilizado por los sismólogos para indicar la atenuación de la respuesta. 


Para índices de amortiguamiento bajos 




1 

2C 


Q 


(1.3-16) 


Por otra parte, el factor de calidad, Q, puede relacionarse con la Capacidad Específica de 
Amortiguamiento, definida en el apartado 1. 2.2.3, como la relación entre la energía 
disipada en cada ciclo del movimiento y la energía almacenada en el sistema; en efecto, 
teniendo en cuenta la relación (1.2-28), y para índices de amortiguamiento bajos puede 
escribirse 

= — (1.3-17) 

E Q 


Las relaciones de frecuencia correspondientes a los puntos R) y R2 para los que el factor 
de amplificación es igual a Q/V2, reciben el nombre de puntos de semipotencia y la 
diferencia entre sus frecuencias asociadas, ancho de banda de sistema. 



Fig.1.29 
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Los puntos de semipotencia, Ri y R 2 , pueden obtenerse al imponer la condición 

1 


es decir, 


8 st V2 2V2C 

1 1 


2V2 C 7(!- P 2 ) 2 +4^ 2 p r 

ecuación cuya resolución proporciona para p 2 los siguientes valores 

P 2 =(1 - 2^ 2 ) ± 2C,y[l + C? 

Rf =(l-2^ 2 )-2cVl + C 2 
r 2 2 =(i-2C) + 2(;^T7c 


por lo que 


y en consecuencia 


si se toma 


el ancho de banda será 


r!-r, ! =4í;Vi + ¡; ! = 


<0 2 - O)? 

O) 2 


0 ) = 


© 2 + ©J 


©, 


-©, =<ü2(;Jl + C 
y para índices de amortiguamiento bajos 

© 2 - ©, = © 2 

siendo entonces el Factor de Calidad 


Q = 


1 


© 


2 C, © 2 - ©, 


(1.3-18) 


(1.3-19) 


(1.3-20) 


(1.3-21) 


(1.3-22) 


lo que sugiere, de nuevo, un procedimiento para estimar el índice de amortiguamiento 
viscoso de un sistema tal y como se indica en el ejemplo E.l-10 
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1.3.3.- FUERZAS TRANSMITIDAS Y AISLAMIENTO DE VIBRACIONES 

Cuando, por ejemplo, un motor está situado sobre una base, transmite a la cimentación o 
estructura soporte fuerzas procedentes de desequilibrios internos o excitaciones externas 
transmitidas por otros sistemas mecánicos. 

Esta fuerza transmitida puede reducirse mediante aisladores con un diseño apropiado de 
sus características. 

Un modelo simple consiste en una base fija (B) sobre la que descansa el motor o 
máquina de masa me, unida a la cimentación mediante un muelle de constante k y un 
amortiguador c. El motor se considera sometido a una fuerza de excitación 

F(t) = F„ sen © 0 1 


1 F (*> 



Fig. 1.30 

en estas condiciones, el comportamiento dinámico del sistema responde a la ecuación 

my + cy + ky = F 0 senco 0 t (1.3-23) 

siendo la respuesta del sistema del tipo 

y = Y sen (® 0 t-<j)) 

La fuerza transmitida a la cimentación será la suma de 
• La transmitida por el resorte 

F k = k Y sen (ü> 0 t - <j)) (1.3-24) 
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® La transmitida por el amortiguador 

F c = c co 0 Y sen (co 0 t - <{>) 
por tanto, la Fuerza Transmitida será 

F k + F c = k Y sen (co 0 t - (J>) + cco 0 Y eos (co 0 t - <])) 
y la amplitud, F t , de la misma 


recordando que 


puede escribirse 


por otra parte 


por lo que 


F t =V(kY) 2 +(cG) 0 Y) 2 
- Y-^/k 2 + (cco 0 ) 2 



F t =kYJl + 


c 2 co 2 


= k YJl + 


4£ 2 km co 2 


co 


o \ 2 


= kY^l + ( 24 ^) 


Y = 


F 0 /k 


®0 s 2 , 4 Í -2 ®0 


í(l - l) +4^ 


(D 


0) 


1 + 4C 


F.=F,- 


2 ®0 
o 2 


CO 


(1 r ) 2 + 4 C 


co 


2 ®n 

co 2 


(1.3-25) 


(1.3-26) 


(1.3-27) 


(1.3-28) 


(1.3-29) 
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y en consecuencia el coeficiente de transmisibilidad, razón de amplitudes entre la fuerza 
transmitida y la fuerza de excitación será 


1 + 4 C 


2 


co 


( 1 - 




(1.3-30) 




r + 


0 + 


la relación de frecuencias, p , para la cual T R = 1, se obtendrá al imponer la citada 

condición; es decir 


1 + 4^ 2 P 2 =(l-p 2 ) 2 +4^ 2 p 2 


por lo que 


P = V2 

T R =1=> (1.3-31) 

p=o 


En la figura 1.31, se recoge, para distintos índices de amortiguamiento, la representación 
gráfica del coeficiente de transmisibilidad correspondiente a la ecuación 1.3-30. 

Como puede observarse, en la zona de frecuencias próximas a la de resonancia, el 
coeficiente de transmisibilidad, Tr, y en consecuencia la fuerza transmitida, son muy 
grandes, sobre todo para índices de amortiguamiento pequeños. 

Por otra parte; la presencia de amortiguamiento disminuye, siempre, la amplitud del 
desplazamiento, pero sólo disminuye la amplitud de la fuerza transmitida cuando la 
relación de frecuencias p < \/2 . Por el contrario, si p > f2 , la presencia de 
amortiguamiento, manteniendo T R < 1, siempre aumenta la amplitud de la fuerza 
trasmitida. 
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COEFICIENTE DE TRANSMISIBILIDAD 



Fig.1.31 


En consecuencia, si la frecuencia de excitación varía, los amortiguadores deberán 
diseñarse de manera que su índice de amortiguamiento sea 

• Lo suficientemente grande como para limitar la amplitud del 
desplazamiento y fuerza transmitida cerca de la resonancia. 

• Lo suficientemente pequeño como para no aumentar la fuerza transmitida 
cuando (1 > /2 
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1.3.4.- RESPUESTA A UNA CARGA ARMONICA DE TIPO GENERAL 
En este caso la fuerza de excitación aplicada sobre el sistema es de la forma 

F(t) = P 0 e iat 

siendo en consecuencia la ecuación del movimiento 


y + 2a^coy + co 2 y = — e 

m 


o „iat 


que tendrá una solución particular del tipo 


y p (t) = Ye' 


iat 


Sustituyendo 1.3-34 en 1.3-33, se obtiene para Y 


Y = 


O 

m 


o bien 


Y = - 


(ca 2 - a 2 ) + 2^a<D¡ 

ni 


co 


(1 — ~j) + 2 ^ ~ ¡ 

co co 


2 k 


y teniendo en cuenta que co = 


m 


Y =- 


k 


1(1 - ~~ r ) 2 + ( 2 C “ ') 2 


co 


co 


es decir Y, es un número complejo que puede expresarse en la forma 


Y = - 


k 


(l-^-) 2 +(2C-i) 3 


„ a . . a 
(1 r ) + «2C- 

CD CO 


CO 


CO 


(1.3-32) 


(1.3-33) 


(1.3-34) 


(1.3-35) 


(1.3-36) 


(1.3-37) 
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o bien 


Y= Ye 


siendo Y I el módulo 


Y =- 


¡ 4> 

k 


(1-^V) 2 +(2C- i) 2 


y <|> el argumento, de manera que 


co 


a 


(ú 


tag* = ® 


1 


a 


o* 


por lo que la solución particular será 


y p (t) = |Y|e i(0 *-* ) 

La ecuación 1.3-36, también puede escribirse en la forma 



1 


a 


(1 r ) + 2C 


0) 


a . 
co 


- H(ia) 


1.3.5.- RESPUESTA A UN IMPULSO 


En este apartado se realizará el estudio de la respuesta de un sistema 
solicitación cuyas características se indican en la figura 1.32 



(1.3-38) 


(1.3-39) 


(1.3-40) 


(1.3-41) 


(1.3-42) 


frente a una 
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FASE I : t < ti 

no hay movimiento, por lo que la respuesta será 

y(t) = o 


FASE II : Ti < t < T 2 Vibraciones forzadas 

La ecuación de comportamiento dinámico del sistema será 

my + cy + ky = F 0 

Cuya solución, tal y como se indico en el apartado 1 .3.1, es de la forma 

y(t) = — + e“ Cffl( '“ Tl) [A eos <» d (t - t, ) + P senco d (t — T t ] 
k 

e imponiendo como condiciones iniciales del movimiento para esta fase 

y(t = Xl ) = 0 
y(t = T,) = 0 


se obtiene 


F 

y(t)=~r 

k 


1 - e 


-CtD(í-X,) 


eos m d (t - t, ) + . ■ senco d (t - t, 


(1.3.43) 


ecuación que para t = T 2 proporciona 


* y 2 =y(t = *2) 

F, 


= “ [l ~ e C(Tl_t,) [cOS © d (T 2 -T,) 


(1.3-44) 


+ 




rsen co d 


( t 2 - x J 


* 


v 2 =y(t = x 2 ) 
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FASE III : t > T 2 Vibraciones libres 

En este caso, la ecuación de equilibrio dinámico será 

my+cy+ky=0 

cuya solución es 


y(t)=e 00(1 Ti) [Acos co d (t - t 2 ) + P sen co d (t-x 2 )] 

realizando el cálculo de las constantes A y B, teniendo en cuenta que las condiciones 
iniciales para esta fase son 


por lo que 


y(t = x 2 ) = y 2 
y(t = x 2 ) = v 2 


y(t)=e~ Cffl( ‘ * [y 2 eos o¡> d (t-x 2 ) + 


(1.3-45) 


A partir de lo anterior, es posible obtener la respuesta a un impulso de breve duración; es 
decir, cuando: xj — > x 2 — » t ; por lo que 


Ar = r 2 - r, -» 0 


con la condición previa, expresada en (1.3-44), se deduce que y 2 = 0, y por tanto, la 
respuesta a un impulso será 


dy(t) = e ía ’ (t T) — -sen co d (t — x) (1.3-46) 

® d 

obtenida al hacer y 2 = 0 en (1.3-45), y donde, dv 2 , representa el cambio brusco de velocidad 
producido como consecuencia del impulso F(x) dx, es decir 

dv 2 = “F(x)dx 
in 


siendo, finalmente, la respuesta 


dy(t) = e ?0>(( T) sen o) d (t-x) 

mco d 


(1.3-47) 
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1.3.6.- RESPUESTA A UNA SOLICITACION DINAMICA DE TIPO GENERAL 


Considérese un sistema ante una solicitación como la que se indica en la 

figura 1.33 



La respuesta en un instante t cualquiera, puede obtenerse por superposición d ellos 
obtenidos a impulsos de duración dx, en los instantes precedentes a t; en efecto, teniendo en 
cuenta lo obtenido en el apartado anterior, puede escribirse 


y(t) = 


i t 

f e -? m( *- i:) F(x)sen co d (t-x)dx 

m ®d o 


(1.3-48) 


la expresión anterior es denominada segunda Integral de Duhamel y representa la 
respuesta de un sistema a una solicitación de tipo general F(x). 

Para sistemas no amortiguados, C, = 0, por lo que 


y(t) = — 

meo 


t 

|F(x)sen ro(t - x)dx 
0 


(1.3-49) 
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1.3.7.- RESPUESTA A SOLICITACIONES SISMICAS 

El diseño de estructuras resistentes a sismos, requiere el conocimiento del movimiento 
producido en el suelo durante el terremoto; dicho movimiento consiste en traslaciones en 
cualquier dirección, combinadas con rotaciones alrededor de cualquier eje. 

Es usual en ingeniería definir la carga sísmica mediante un acelerograma, que 
representa las componentes, en las direcciones horizontales y verticales de la aceleración 
producida, en un determinado lugar de la superficie del suelo, por un terremoto. 

En la figura 1 .34, se representa la componente N-S de la aceleración del suelo obtenida 
durante el terremoto de El Centro (1940); dicha aceleración puede integrarse obteniendo de 
esta forma los valores de la velocidad y del desplazamiento del suelo. 


0.20 g. 
0.10 


ACELERACION DEL SUELO 



0.10 

0.20 


L 34 cm/seg 


VELOCIDAD DEL SUELO 



TIEMPO (seg) 


Fig.1.34 

En ocasiones se precisa disponer de acelerogramas adaptados a las características del 
emplazamiento, no existiendo acelerogramas reales ad-hoc, por lo que se recurre a generar 
registros artificiales. En la figura 1.35, se representa el acelerograma sintético de Bogdanoff 
(1968), que se utiliza en ocasiones para generar, en sistemas lineales, desplazamientos 
semejantes a los reales; utilizando para ello una formula de representación de la aceleración 
del suelo dada por la expresión. 


10 

a g (t) = te' 0333, ^cos(co j t + (p j ) ;t>0 


j=n 
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los valores de coj y (pj ,se indican en la tabla 1.3.1 


j 

(Dj 

<Pj 

1 

6,00 

3,7663 

2 

8,00 

1,3422 

3 

10,00 

4,8253 

4 

11,15 

0,2528 

5 

12,30 

4,5204 

6 

13,25 

1,8834 

7 

14,15 

1,3320 

8 

16,20 

1,7852 

9 

17,35 

0,1517 

10 

19,15 

2,4881 


Tabla 1.3.1 



En el apartado 1.4.6, se vio que la respuesta forzada de un sistema con 1 g.d.l. a una 
carga F(t) venía dada por la expresión 


! t 

x(t) = fe^ <D(, ” T) F(x)sen co d )(t-T)dT 

meo, J 


(1.3-50) 


J d 0 
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En el caso de un desplazamiento x s (t) y otro relativo x(t), respondiendo el 
comportamiento a la ecuación 

m[ a s (t) + x(t)] + c x(t) + k x(t) = 0 

o lo que es igual 

m x(t) + c x(t) + k x(t) = - m a s (t) (1 .3-5 1 ) 


# x s<*> tK— X(t> —7^ 



* X s<*) tP 


Fig.1.36 

Por lo que la carga equivalente a la aceleración del suelo a s (t) será 

F(t)=-ma s (t) (1.3-52) 

sustituyendo la relación (1.3-52) en (1.3-50) se obtiene que la respuesta del sistema 

1 t 

x(t) = íe“ c< “ (tt) a s (T:)sen© d (t-x)dT (1.3-53) 

®d o 


como puede observarse la cantidad subintegral depende del terremoto dado a través de la 
aceleración del suelo a s (t), y de la estructura escogida mediante la frecuencia propia © y el 
índice de amortiguamiento C,. 
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1.3.8.- TRATAMIENTO NUMERICO 


En la mayoría de los casos, y en especial con la solicitación sísmica, las cargas actuantes 
tienen una definición complicada que impide la obtención de soluciones de tipo analítico. 
Por ello se precisa disponer de un método de solución numérica que, de forma 
suficientemente aproximada, garantice la obtención de la respuesta. 

En lo que sigue se va a representar un método mixto en el que se utiliza una solución 
exacta con una aproximación de la ley real. Es decir, se aproxima el acelerograma mediante 
trozos de recta y se integra exactamente el acelerograma modificado. El error de 
procedimiento se concentra así en la fase de modelado de la solicitación que, por otro lado, 
suele almacenarse de esta forma en los observatorios sismológicos. 

El desarrollo se realiza de manera que las formulas obtenidas sean fácilmente 
programables. 

Supóngase que la aceleración del suelo, x s (t) = a g (t), puede representarse mediante 
una función lineal tal y como se indica en la figura 1.37 



Es decir 

a g (t) = a ¡ + 7^ : t¡ <t<t¡ +1 


A a ¡ = a ¡ + i “ a ¡ 
At¡ = t ¡+1 - t¡ 


o bien designando por 
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puede escribirse 


a g (t) = a¡ +^ L ( t - t ¡) 


(1.3.54) 


por lo que la ecuación diferencial (1.3-51) toma la forma 

Ao 

x(t) + lífo x(t) + co 2 x(t) = -a¡ - — d- (t - t. ) 

At, 


(1.3-55) 


La solución de la ecuación diferencial (1.3-55), será la suma de la solución de la 
homogénea 

x h (t) = * [A eos ® d (t - t¡ ) + B sen co d (t - 1¡ )] 

más una particular de la completa; que será, para este caso, una expresión lineal 

x p (t) = at + p 

por lo que 

x p (t) = <x 
* p (t) = 0 

valores que sustituidos en (1.3-55) proporciona 

2^íoa + o) 2 (at + P) = -(a i - 

1 At, 1 At, 


o bien 


(ff> 2 p + 2^oaa) + co 2 at = -(a, -^^-t.) --^-t 

At, 1 At, 


identificando coeficientes, se obtienen los valores de a y p 

Aa, 


* <n 2 a = — - 


At, 


1 

Aa¡ 

^ 2 
co 

At, 
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<0 


, Aa- 

co P + 2C<oa = -a¡ H 1, 

' At, ' 

Aa¡ 2C Aa, 

-a. + — +— — — 
At, (¡o At, 



J_Aa, K 

CO 2 At. (D 


Siendo entonces la solución de la particular 
x p (t) = cxt + p 

o 2 At¡ co co At, co 
co 2 co 2 At, 1 co 


La solución a la ecuación (1.3-55), para t¡ <t<t¡ +1 , es por tanto 

x(t) =e ■^®(‘-‘i)[Acos(o ( ,(t-t i )+Bsen(o d (t-t i )]- 

co 2 co 2 At, 1 co 


siendo co d =(Oyjl-C, 2 . 

Las constantes A y B se calcularán imponiendo las condiciones: 


t = t, => 


X = X, 


X = x. 


siendo el desplazamiento, para t = t¡ 


x (t = t,) = A l T + 


2? Aa, 


= x, 


co 2 co 3 At, 

expresión, a partir de la cual es posible determinar el valor de A 


a . Ii+ ±*í a *. 


co 2 co 3 At¡ 


(1.3-56) 
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por otra parte, al ser la velocidad igual a 

x(0 = _ ^®e ' c<o ^ [A eos © d (t - t¡ )+ Bsen co d (t - t¡ )]+ 

+ e -«•(•-*■) [_ A<» d seno d (t - 1¡ )+B(o d cosco d (t - t, )] — (1.3-56b) 

1 Aa¡ 
o 2 At, 


el valor que toma para t = t. será 

x(t = t¡) = -C;a>A + G) d B- 
de donde puede obtenerse el valor de B 


1 Aa, 


J co 2 At, 


• = x, 





B= — 

. . 2^-lAa, C 

x¡ +C® X ¡ ^ ~+—a¡ 


®d 

© 2 At. co 



En el supuesto de un intervalo de tiempo constante, At. = cte = At , haciendo t = t i+1 , 
se obtienen para A y B los siguientes valores 


. a¡ 2C a... -a, 

A = x. + — L = x i +a¡ 

co 2 co 3 At 11 


1 + - 2í ; 


(O 2 (0 3 At 


2C, a 


i+l 


co At 


B = — 

©,, 


©, 


. „ 2^ 2 -1 / v £ 

x í+C® x í - a J + T a ¡ 

x¡ +^©x¡ 


© At 

2 ¿; 2 -i 


2 a a a i+l + a ¡ 


© At 


2 C 2 -i + C 


y para la solución particular 


a¡ 1 Aa¡ 


© 2 © 2 At. 

a, 


t-t, - 


1 a ¡ + i- a ¡ 


© 2 © 2 At 

2 C 

2T7 a i +a ¡ + i 
© At 


At 


L 

! a t 

s» 


© J 


_ 2 C 

i= 


© 


zl + JL 

© 2 © 3 At 


x P (t) = 
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La sustitución de A y B en la ecuación (1.3-56), proporciona una ley de recurrencia que 
permite determinar el valor del desplazamiento en el instante i +1 


_ -CcoAt 
X i+1 e 


x. +a. 


" 1 

1 

2C a ,+, 

(O 2 

co 3 At 

co 3 At 


eos co d At + 


+ e 


-C <D At 


(Oh 


2C 2 -1 

x ¡ +C®x¡ — a i+ i -a, 


oo At 


2C 2 -1 , C 

co 2 At co 


senco d At- 


2C, 

a ¡ 3 . + a ¡+i 

co At 

y reagrupando términos 


2 S 1 


<o 3 At co 2 


(1.3-57) 


i.+.= e' C ’* J 


cosco d At + 




rsen co. At lx¡ + ( e ® At 


sen co. At 


(Oh 


•x. + 


+ (e 


co At 


í V f 1 

eos co d At + 

^co co 3 At; 

V 


2C -1 , C 

co 2 At co 


sen co d At 


®h 


2£ 

<o 3 At I 


a, + 


+ (e 


co At 


2C . 2C 2 -1 senco d At 

-coscOh At + 


co At 


co At co H 


2C 1 


•a., 


1 i i 2 **i+l 

co At <o 


(1.3-58) 


finalmente, la velocidad en el instante i +1, se obtendrá al sustituir las constantes A y B en 
la expresión (1.3-56b) 


-C ©At 


CO 




rsencOn Atlx, +e " ímAt 


+ < 


-C, (o At 


O 

— 2 = — +— 

v co At co y 


cosco d At 


eos co d At- 


sen co d At 




:senco d At 


+ 


2 « + 1 


co 3 At co 2 


(co d sen <o d At - ^co co d At) 


<o At 


ca. + 


+ <e 


-C, co At 


2 - f 


2^-1 

co 2 At 


cosco d At- 




- sen co d At 


(co d sen co d At - Cfü co d At) 

co At 


+ - 


co 2 At 


?-a 


i+l 


(1.3-59) 
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Las ecuaciones (1.3-58) y (1.3-59), en forma matricial quedan como sigue 


x ¡ + l" 

= D 

x ¡ 

+ E 

a ¡ 

>¡ +l . 


> ¡. 


_ a ¡+i. 


y si las matrices D y E se escriben en la forma 


D = 

d.. 

d,2 

• 

E = 

e n 

e i2 _ 

— 

_ d 2, 

d 22. 


— 

_ e 21 

e 22. 


los coeficientes de las mismas serán 
❖ Para la matriz D 


d u = e <mM 


coso. At + -Jg^ rsen cd, At 


^ _ „ -C »A. Sen ®d 

U 12 _ e 


CO,, 


co 


d 21 =-e-^ At 




senoo d At 


d 22 =e \ 


eos co d At - 


c 


a / i 1 ? 


rsenoo,, At 


♦ Para la matriz E 


-C co At 


f 1 

— + — — eos oo ri At + 

co 3 At J d 


^2 C, 1 -1 £") sen oo d At 


o At (o j 


2£ 

co 3 At 


e i2 = e “ At 


2C . 2C 2 -1 sen oo d At 

^ costa. At + ^ d 


- . 

o At 


co At oo rf 


+ 


2C 1 


oo 3 At © 2 


e 2i =e 


-C, co At 


r 2C~l C 
— +— 
v 00 At 00 J 


cosoo d At - 


sen ro d At 


( 2^ 1 

^oo 3 At oo 2 


(oo d sen CD d At - t^oo co d At) 


+ 


oo 2 At 


© At 


2C 2 -1 
oo 2 At 


cosoo d At- 




sen oo d At 


— (oo d sen co d At - £co oo d At) 


G) At 


+ - 


1 


or At 
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Si el índice de amortiguamiento = 0, entonces ío d = o>, por lo que los coeficientes de 
las matrices D y E se simplifican, tomando los siguientes valores 

♦ Matriz D 



4 Matriz E 
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1.3.9.- INSTRUMENTOS DE MEDIDA DE LAS VIBRACIONES 


Son aplicaciones particulares de sistemas con 1 g.d.l., y consisten fundamentalmente en 
un sistema masa-muelle-amortiguador, montado sobre una base soporte tal y como se 
indica en la figura 1.38 



Fig. 1.38 

Si se designa por x (t) al desplazamiento de la masa m ; el desplazamiento relativo u(t) 
entre la masa y la base será 

u(t) = x(t) + u s (t) (1.3-60) 

siendo la ecuación que caracteriza el movimiento del sistema 

mü(t) + cx(t) + ku(t) = 0 (1.3-61) 

o lo que es igual, teniendo en cuenta (1.3-60) 

m ü(t) + c x(t) + k u(t) = -mii s (t) (1.3-62) 

y dividiendo por la masa 

ü(t) + 2£oo ú(t) + o) 2 u(t) = -ü s (t) (1.3-63) 


suponiendo que el movimiento del soporte, u s (t) , es armónico de frecuencia, co 0 , es decir 

u s (t) = U s senco 0 t (1.3-64) 
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se tiene 


ü(t) + 2^co ú(t) + oo 2 u(t) = -U s G) 2 sen oo 0 t (1.3-65) 


siendo la solución a la ecuación anterior 

u(t) = U sen (cd 0 t - <j>) 


(1.3-66) 


donde la amplitud, U , y el ángulo de fase, <¡>, pueden expresarse mediante las relaciones 


2£ — 
<j) =arctan ^r- 

(D 2 


(1.3-67) 


U 


uA 

CO 


(l_°^) 2 +4 ^^ 

CO £0 


(1.3-68) 


por otra parte, teniendo en cuenta que 

ü,d)= ",(0 ^ Ü,= co. ! U, 

la relación de amplitudes entre el movimiento relativo y el del soporte, puede escribirse en 
la forma 


IJ 

Ü s 



co 


ca 


~) 2 +4C 


co 2 


(1.3-69) 


Las expresiones 1.3-68 y 1.3-69, constituyen la base del funcionamiento de los 
instrumentos denominados VIBROMETRO y ACELEROMETRO; en efecto 

♦ Si la frecuencia natural co, del sistema es pequeña, se tendrá que co 0 /co— »<», y por tanto 


lim -^- = 1 => U = U S 

<»o/®->'» u s 

dicho de otro modo, el desplazamiento relativo es igual al desplazamiento del soporte; por 
lo que el registrador de vibraciones lo que realmente está midiendo es el movimiento del 
soporte. 
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Este instrumento de medida recibe el nombre de VIBROMETRO, y las bajas 
frecuencias naturales que lo caracterizan, implica que estos sistemas sean excesivamente 
voluminosos, motivo por el que no suelen utilizarse en la mayoría de las aplicaciones 
prácticas. 

En la figura 1.39 se representa la variación de U/U s y el ángulo de fase <j> en función de 
la relación de frecuencias <ü 0 /(d 



Fig.1.39 

♦ Por el contrario para frecuencias naturales altas a> o /a>-»0, y la ecuación 1.3-69 se 
transforma en 



es decir, el registrador de vibraciones mide la aceleración del soporte, por lo que recibe el 
nombre de ACELEROMETRO. 

En la figura 1.40, se representa, en función de la relación de frecuencias a> 0 /co , la 
variación de la expresión 


1 
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Fig. 1.40 

La figura ha sido dibujada para una distorsión máxima de la medida del acelerómetro del 
6 %, y en ella puede observarse que el límite de utilización del mismo disminuye al 
aumentar el índice de amortiguamiento. 


1.3.10.- ESPECTRO DE RESPUESTA ELÁSTICO 


Aunque el concepto de espectro de respuesta elástico tiene orígenes anteriores (Biot 
1941) y es utilizado en otras ramas de la ingeniería (Espectro de Choque), la comprensión 
de su utilidad y los primeros intentos por utilizarlo, como medio de definición de la acción 
sísmica, fueron debidos a Housner(1941); quien pretendía encontrar una medida de la 
intensidad del terremoto relacionada con la respuesta de las estructuras de forma racional y 
no subjetiva. Por ello propuso estudiar la respuesta máxima, a un terremoto dado, de 
sistemas con un grado de libertad y diferentes características dinámicas; reduciéndose los 
parámetros del sistema, en este caso a dos : 


♦ La frecuencia propia 



que mide la relación entre la masa y la rigidez del sistema 
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♦ El índice de amortiguamiento 

^ _ c _ c _ c 
C crit 2^¡km 2m© 

que relaciona las propiedades de disipación de la energía con aquellas. 

En el apartado 1.3.8, se ha presentado un procedimiento que permite determinar la 
respuesta x (t) de un sistema ante una solicitación sísmica caracterizada mediante un 
acelerograma a (t). 

En la figura 1.41, se presenta la respuesta de tres estructuras sometidas a una excitación 
sísmica, cuyo acelerograma está representado en la parte superior; en todas ellas se ha 
considerado un índice de amortiguamiento C, = 2%, la diferencia entre las tres se refiere a 
su frecuencia natural. 


EL CENTRO 
COMPONENTE SOOE 




PERIODO NATURAL DE VIBRACIÓN 
T (seg) 


Fig. 1.41 
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Sin embargo, desde un punto de vista práctico solo se buscan los valores máximos, 
independientemente del signo, es decir 

X máx = máx | x(t) | 

en la parte inferior de la figura 1.41, se ha dibujado el valor máximo de la respuesta en 
desplazamientos, en función del periodo natural de la estructura. 

El gráfico corresponde a lo que habitualmente se denomina ESPECTRO DE 
RESPUESTA EN DESPLAZAMIENTOS, y se designa por SD. 

SD = máx|x(t)j (1.3-70) 

De forma análoga a lo descrito previamente podría determinarse el ESPECTRO DE 
VELOCIDAD ( SV ) y el de ACELERACIÓN ABSOLUTA ( SA ). 

SV = máx|x(t)| (1.3-71) 

SA = máx| x(t)+a s (t)¡ (1.3-72) 


El ESPECTRO DE RESPUESTA ELÁSTICO es, por tanto, un diagrama de 
representación de los valores máximos de la respuesta de un sistema amortiguado con un 
grado de libertad, sometido a una excitación dinámica, en función de la frecuencia natural 
del sistema o alguna magnitud relacionada con ella. 


Las expresiones analíticas correspondientes al espectro de respuesta son las siguientes: 

♦ Desplazamientos ( SD ) 


SD = 


-| t 

— . fa(x)e“ ítD ” (,_T) sen c o n ^l-C, 2 (t-T)dt 

nVl^O 


(1.3-73a) 


máx 


♦ Velocidad ( SV ) 


SV = 


- |a(x)e ?c ° n(t T) eos (o n ^í-C 3 1 (t-x)dx + 


+ 


C i 


VT-í 2 


ja(T)e e<D " (t T) sen (t - T)dx 


(1.3-73b) 


máx 



1-79 


❖ Aceleración ( SA ) 


SA: 


< ° n ^ - Ja(t)e íc ° n(t T) sen oo n -y/l - C, 2 (t - x)dx + 

Vl-C 2 o 


+ 2«) n ^|a(T)e ?<D " (t T> eos (o n ^l-C, 2 (t-x)dx 


(1.3-73c) 


máx 


Por otra parte, teniendo en cuenta que la integral 

t 

Ja(x)e~ í; ' B '’ (t ~ i;) sen <a n ^l-C ¡ 2 (t - x)di: (1.3-74) 

o 

tiene dimensiones de una velocidad, se designó con el término PSEUDOVELOCIDAD ( 
PSV ), al máximo del valor absoluto de la citada integral, y con objeto de utilizar una 
escala independiente de la frecuencia, se definió el espectro de respuesta en 
desplazamientos, para estructuras con bajo amortiguamiento (co d « co) , como : 


PSV 

SD = — — (1.3-75) 

co 

La magnitud PSV, está relacionada con la energía máxima almacenada en el sistema: 


E 


máx 



[ x (t)]Lx = 
"psv ] 2 

co 

[psv] 2 


(1.3-76) 


Otra medida del desplazamiento máximo se obtiene mediante la relación 


SD = — — (1.3-77) 

co 

La cantidad PSA se denomina ESPECTRO DE RESPUESTA DE 
PSEUDOACELERACIÓN. En las expresiones (1.3-73) puede observarse que sí C, = 0; 
PSA = SA, mientras que, PSV * SV. 
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En la figura 1.42, se recogen, para la componente SOOE del terremoto de El Centro (£, = 
5%), los dos espectros de aceleración referenciados. 



Fig. 1.42 

El PSA tiene unidades de aceleración y está relacionado con la fuerza máxima en el 
muelle 


F -1. = k l x (0L = k SD (1.3-78) 

y según (1.3-77) 

PSA 

^máx - k— p = m PSA (1.3-79) 

GD 


Se obtiene así una fórmula del tipo “masa por aceleración”, que permite definir la 
fuerza estática equivalente al movimiento sísmico, ya que produce el mismo 
desplazamiento máximo. 


Esta idea es la base de los métodos de fuerzas equivalentes contenidos en las Normas y 
Códigos de Diseño. 
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Los tres espectros de respuesta: de desplazamiento ( SD ), de pseudo velocidad (SV ) y 
de pseudo aceleración ( PSA ), están relacionados mediante las expresiones: 


PSV = ro SD 


PSA = ro PSV 


(1.3-80) 


y son formas diferentes de representar la misma información referente a la respuesta del 
sistema. Como cualquiera de los tres espectros puede ser obtenido del otro, es posible 
realizar una representación conjunta en el denominado Diagrama Tripartito, que en escala 
logarítmica permite dibujar las magnitudes espectrales ( SD, PSV, PSA ) en función de la 
frecuencia natural, ro, ó el período, T. 

Tomando logaritmos en (1.3-80) 


log PSV = logSD + logro 


log PSA = log PSA - log ro 


(1.3-81) 


estas relaciones, en el diagrama ( log ro - log PSV ), son rectas paralelas a las bisectrices del 
primer y segundo cuadrante, véase figura 1.43 



Fig. 1.43 



VEL0C1TY (in./sec) 
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En la figura 1.44, se recoge el diagrama tripartito del terremoto californiano KERN 
COUNTY. 


KERN COUNTY, 21 de Julio de 1952 
Magnitud 7.7 Aceleración máxima del suelo 0.057 g 
PASADENA - CALTECH Componente SOOE 
t-0,2, 5, 10 y 20 % 



g. 1-44 


VELOCITY, ¡n./sec 
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En la figura 1.45 se recoge el diagrama tripartito de el terremoto de El Centro, (1940) 



Rcsponsc Spcctra, El Centro Earthquakc, May 18, 1940, 
North-South Dircction, 


Fig. 1.45 
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1.3.11.- ESPECTROS DE DISEÑO 


Sin embargo, la utilidad del espectro de respuesta no se limita a la caracterización de un 
terremoto individual. 

Conocidos los espectros para diferentes terremotos y a la vista de la forma semejante 
que presentan, se intentó buscar una curva promedio. Para conseguirlo, fue necesaria la 
normalización de los espectros, lo que se llevó a cabo mediante dos procedimientos 
diferentes. 

El primero , históricamente, se debe a Housner (1959), que observó que para una 
banda ancha de frecuencias el espectro de velocidad era horizontal, por lo que propuso 
utilizar como base de normalización la denominada Intensidad Espectral 

IS = JJ = 0 2 ’ 5 PSV( co) doo (1 .3-82) 

o área comprendida bajo el espectro de Pseudo-Velocidad entre los periodos 0,4 y 2,5 
segundos; recuperando el espectro original multiplicando por el valor de la intensidad 
espectral S 0 ; en la figura 1.46 se muestran los primeros espectros utilizados para el 
proyecto sísmico de estructuras. 

Housner utilizó como base de la normalización los registros correspondientes a los 
terremotos de El Centro (1934), El Centro (1940), Olimpia (1949) y Taft (1952). Los 
valores de S 0 ,se indican en la siguiente tabla 


Earthquake 

Direction 

So 

Peak ground 
acceleration 

El Centro 

E.W 

2.2 

0.08 

1934 

N.S 

2.2 

0.26 

El Centro 

E.W 

2.8 

0.22 

1940 

N.S 

3 

0.33 

Olympia 

E.W 

2.1 

0.31 

1949 

N.S 

1.8 

0.19 

Taft 

E.W 

1.8 

0.16 

1952 

N.S 

1.7 

0.18 


Tabla 1.3.11-1 





PeriocJ (sec) 

Housner acceleration design spectra, peak ground acceleration = 0,125 g 

Fig. 1.46 
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La búsqueda de un factor de normalización más sencillo condujo a la utilización de la 
aceleración máxima del suelo, a máx , por lo que la expresión 1.3-73a puede escribirse como: 


X 


máx 


máx 


co 


í-^^-e í<0 " ( * T) sen co n -^/l-^ 2 (t-T)dx 

l nA /l-^ 2 0 «máx 


(1.3-83) 


máx 


Un espectro típico, puede verse en la figura 1.47 que corresponde de nuevo al terremoto de 
El Centro, y en él puede observarse que para T = 0 el PSA arranca del valor unidad, ello se 
comprende teniendo en cuenta que periodos nulos se corresponde con estructuras 
infinitamente rígidas por lo que la aceleración es la del suelo. 



Acceieratíon, veiocity, ana dispiacement ampilfica- 
tions plctted as a function of oerlod for 5% damoing for the SOQE 
component of El Centro, ¡ne Imperial Vailey eartftquake of May 

18 . mo. 


Fig.1.47 



1-87 


Cuando T— > oo, se produce una situación semejante pero con los desplazamientos, 
debido a la absoluta flexibilidad del muelle la masa permanece en su lugar por lo que el 
desplazamiento relativo coincide con el del suelo. 

Las ideas anteriores, extrapoladas para cualquier espectro, condujeron a la segunda 
aproximación que con idéntico objetivo fue realizada por Newmark en una serie de 
artículos publicados entre 1965 y 1973. 

Newmark y sus colaboradores estudiaron los espectros de respuesta para diferentes tipos 
de impulsos comparándolos con los obtenidos para terremotos con objeto de encontrar 
reglas subyacentes. 

En este sentido, observaron que, en un diagrama tripartito, la forma general de un 
espectro promemediado tomaba la forma que se indica en la figura 1.48, donde pueden 
observarse zonas claramente diferenciadas: 




Frequency, f (Log scale) 

Genera) shape of a smoothed response spectrum. 


Fig.1.48 

❖ En el rango de bajas frecuencias, la respuesta máxima en desplazamientos coincide 
con el desplazamiento del suelo. 

❖ Para altas frecuencias, la aceleración máxima es prácticamente igual a la 
aceleración del suelo. 
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❖ En el rango de frecuencias intermedias puede distinguirse una región de respuesta 
en velocidad constante amplificada respecto a la máxima del suelo, que está 
flanqueada por una zona de amplificación de desplazamientos y otra en 
aceleraciones. 

En la figura 1.49 se recoge nuevamente el diagrama tripartito de el terremoto de El 
Centro, donde puede observarse la irregularidad de los espectros de respuesta para un 
movimiento sísmico. 



Response spectra lor 2. 5. and 10% damping for the SOOE cc'.iponent of El Centio. the Imperial 
Vaíley earthquake ol May 18. 1940. together with the peak ground motion. 


Fig. 1.49 

Utilizando diagramas del tipo indicado en la figura 1.48 es posible promediar resultados 
para los terremotos registrados y obtener así una envolvente de los mismos. Newmark 
propuso utilizar líneas rectas para obtener la envolvente en el diagrama tripartito y con Hall 
sugirió los siguientes valores de los factores de amplificación: 
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FACTOR DE AMPLIFICACIÓN 

ÍNDICE DE 

AMORTIGUAMIENTO 

DESPLAZAMIENTO 

VELOCIDAD 

ACELERACIÓN 

0 

2.5 

4 

6.4 

0.5 

2.2 

3.6 

5.8 

1 

2 

3.2 

5.2 

2 

1.8 

2.8 

4.3 

5 

1.4 

1.9 

2.6 

10 

1.1 

1.3 

1.5 


Tabla 1.3.11-2 

Si para 1 g, la velocidad máxima del suelo se supone igual a 48 in/s y el desplazamiento 
36 in, teniendo en cuenta que para el terremoto de El Centro a s = 0,33 g, con las 
consideraciones previas la velocidad y el desplazamiento del suelo serán: 

v s =48 x 0,33 = 15,84 in/s = 16 in/s 
d s =36 x 0,33 = 11,88 in =12 in/s (1.3-84) 

Por lo que los valores amplificados de la aceleración, la velocidad y el desplazamiento, 
se obtendrán a partir de los factores de amplificación de la tabla. Así, para un índice de 
amortiguamiento de un 2% : f d = 1,8 ; f v =2,8 ; f a =4,3, por lo que para 

d 0 =36in ; v 0 =48in/s ; a 0 = lg , se obtiene: 
a = f a a s = 4,3 x 0,33 = 1,42 g * l,4g 

V = f v v s =2,8 x 15,84 = 44,35 in/s as 45 in/s (1.3-85) 

d = f d d s =1,8 x 11,88 = 21,38 in «22 in/s 


con ello puede representarse el diagrama tripartito indicado en la figura 1.49 
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Frequency, cps 

SMOOTHED TRIPARTITE LOGAR ITH MIC (0,33 g ; 2% CRITICAL DAMPING) 


Fig.1.50 

Tras una serie de estudios estadísticos la Nuclear Regulatory Conmission de USA 
propuso los espectros recogido en la figura 1.51 que forman la base de todos los espectros 
posteriores propuestos por las Códigos de buena práctica. Las frecuencias de control son: 

33, 9, 2.5 y 0.25 Hz y se supone que para una aceleración de 1 g, la velocidad máxima del 
suelo es igual a 48 in/s y el desplazamiento 36 in. 

La tabla recoge los valores de amplificación correspondiente a cada frecuencia, así como 
la ley que permite interpolar la curva correspondiente a cada amortiguamiento 

Las ordenadas corresponden a los valores resultantes de tomar la media más una 
desviación típica, lo que implica que la probabilidad de superación es sólo de un 15,9%. 
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1.3.12.- ESPECTROS DE UTILIZADOS EN LOS CÓDIGOS DE PRÁCTICA 


Si se transforma la figura 1.48 en un diagrama (PSA,T) y se tiene en cuenta que el punto 
de anclaje de la aceleración del suelo se estima en 30 Hz y la transición en 6 Hz, se obtiene 
la figura 1.52 



Fig. 1.52 

Con referencia a la figura 1.48, la zona AB al ser de PSA constante se transforma en 


A B . 

27C 

En CB el PSV es constante y por tanto PSA = coPSV =— PSV, en consecuencia la 
rama B*C* cae según 1/T. 


471 ^ 

Finalmente el tramo CD tiene un SD constante: PSA = co 2 SD = — ¿-SDy por tanto la 
rama C*D* cae según 1/T 2 . 
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La figura 1.53 recoge la forma que presenta el espectro de Newmark al expresar el PSA 
en función de T. 



AEC rcgulaiorv guitic spcctrum (S4 ncrcennlel. 


Fig. 1 .53 


En líneas generales los espectros promedios propuestos por los diferentes códigos 
coinciden con los indicados previamente: se establecen unos puntos de anclaje en períodos, 
se definen unos factores de amplificación y una ley de caída para periodos altos; todo ello 
debido a la relación 1.3-79 y a la utilización generalizada de un método de fuerzas 
equivalentes. 
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El eurocódigo 8 (1988), por ejemplo, establece el espectro elástico normalizado p e (T), 
figura 1.53, para un índice de amortiguamiento £ = 5%, mediante cuatro parámetros: Po, Ti, 
T 2 y k 



Fig. 1.54 

Y las siguientes expresiones: 

0 < T < T, ; Pc ( T ) = l + ^(p 0 -l) 

T,<T<T 2 ; Pe (T )= P 0 (1.3-85) 

ÍT V 

T> Tj ; P«(t)= -i- p, 

V 1 

advirtiendo que esta representación no es válida para valores de T > 5 segundos. 

Los parámetros deben ajustarse según el tipo de terreno y en el caso de suelo rocoso 
toman los siguientes valores: p 0 = 2,5 ; Ti = 0,2 s ; T 2 = 0,49 s y k = 1. 

Para otros valores del índice de amortiguamiento las ordenadas se corrigen con el factor 
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Como puede observarse la forma corresponde a la propuesta por Newmark y Hall, 
suplementada por la transición lineal entre T = o y Tj. 

También el Código Modelo del CEB para Proyecto Sismorresistente de Estructuras de 
Hormigón Armado (1985), propone un espectro como el del EC8 con valores 
recomendados de Ti = 0,2 s ; T 2 = 0,40 s . 

La (NCSR-02) : Norma de construcción sismorresistente: parte general y edificación 
(RD 997/2002, de 2-'Sep) establece un espectro normalizado de respuesta elástica en la 
superficie libre del terreno, figura 1.54, para aceleraciones horizontales, con un índice de 
amortiguamiento del 5%, definido por los siguientes valores: 


t<t a 

; a (T ) = 1 + 1,5 — — 

1 A 


t a <t<t b 

; a(T)= 2,5 

(1.3-85) 

t>t b 

i »(T)= k| 



donde: 


❖ a (T)= Valor del espectro normalizado de respuesta elástica 

❖ T = Período propio del oscilador en segundos 

❖ K = Coeficiente de contribución, que tiene en cuenta la influencia de los 

distintos tipos de terremotos esperados en la peligrosidad sísmica 
del emplazamiento 



TIPO DE TERRENO 

C 

I 

Roca compacta, suelo cementado 0 granular muy 
denso 

1.0 

II 

Roca muy fracturada, suelos granulares densos 0 
cohesivos duros 

1.3 

III 

Suelo granular de compacidad media, 0 cohesivo de 
consistencia firme a muy firme 

1.6 

IV 

Suelo granular suelto, 0 suelo cohesivo blando 

2.0 


❖ C = Coeficiente que tiene en cuenta las características geotécnicas del terreno 
de cimentación 
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❖ T A y T B , períodos característicos del espectro de respuesta de valores: 


T a 

Tb 



2,5 


a (T) - ' 

3.0- 



1 » 

T b periodo de oscilación, t 

Fig. 1.55 



para algunas estructuras y modos de vibración con períodos superiores a 4 segundos, 
siempre y cuando se justifique, se podrán utilizar valores espectrales menores. 
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2.1. -INTRODUCCIÓN 


En el capítulo anterior, se han desarrollado los principales conceptos utilizados 
en el análisis dinámico de sistemas estructurales. Aunque la modelización de un 
sistema mediante un grado de libertad proporciona, en muchos casos, una buena 
aproximación de cara a un análisis preliminar, no siempre aporta resultados 
suficientemente precisos para un diseño final ya que la mayoría de los sistemas 
estructurales tienen particularidades que exigen un modelo mas detallado, lo que 
induce la introducción de nuevos g.d.l.. 

Como ya se mencionó previamente, el número de grados de libertad de un 
sistema es el número de coordenadas independientes que es preciso especificar para 
definir el movimiento de un sistema. 

En un sistema con N g.d.l. cada una de las masas puede moverse 
independientemente de las otras, y solo bajo ciertas condiciones todas ellas tendrán 
un movimiento armónico con la misma frecuencia. Esta frecuencia recibe el nombre 
de FRECUENCIA NATURAL y la configuración deformada del sistema para esta 
frecuencia MODO DE VIBRACIÓN , siendo el número de frecuencias naturales y 
modos de vibración de un sistema igual al número de grados de libertad del mismo. 
Es conveniente, por tanto, pensar en el Modo de Vibración como una configuración 
deformada en la que el movimiento de cada masa es una vibración armónica, 
alrededor de la posición de equilibrio, con una frecuencia especifica asociada a este 
modo de vibración. 


2.2.- ECUACIONES QUE CARACTERIZAN EL MOVTMIU.NTO 


Un sistema con dos grados de libertad solo requieren dos coordenadas para 
describir su movimiento. Así, el indicado en la Figura 2.1., tiene 2 g.d.l. que 
corresponden a los dos movimientos de traslación de la masa m 



Fig. 2.1 
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En la Figura 2.2., se representa otro sistema con 2 g.d.l.; el movimiento se 
describe a partir de las coordenadas x,(t) y x 2 (t), que definen las posiciones de las 
masas nq y m 2 en cualquier instante de tiempo t. 



Fig. 2.2 

La aplicación de la 2 a ley de Newton permite escribir 

* MASA m, 

/Mj i'j = Fj + (x 2 - Xj) + c 2 (x 2 - ij) - k t x 1 - Cj Xj 

* MASA ni 2 

^2 X 2 ~ ^2 ~ ^“3 X 2 ~ ^3 X 2 ~ ^*2 C ^2 — — ^2 C ^2 ~ ^1^ 

es decir 

m l x'j + (Cj + c 2 ) x 1 - c 2 x 2 + (k t + kj Xj - x 2 = F í (2.2-1) 

m 2 *2 ~ C 2 *1 + ( c 2 + ^ 3 ) ^2 ~ ^2 x i + (1% + k 3 ) x 2 = F 2 (2.2-2) 

las ecuaciones (2.2-1) y (2.2-2) forman un sistema de ecuaciones diferenciales 
acopladas, ya que los dos movimientos x 2 (t) y x 2 (t) aparecen en ambas ecuaciones. 
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El sistema de ecuaciones puede escribirse en forma matricial como sigue 


m x(t) + c x(t) + k x(t) = F(t) 


(2.2-3) 


Obsérvese que la ecuación es totalmente semejante a la que caracteriza el 
comportamiento de un g.d.L, salvo por el carácter matricial de sus términos; 

m » c y k , son las matrices de MASA, AMORTIGUAMIENTO Y RIGIDEZ del 

sistema 


m t 

0 

; c = 

C l +C 2 C 2 

; k = 

k í +k 2 ^2 

0 

m 2 


~ C 2 C 2 +C 3 


+ 

i 


y, x(t) , x(t) , x(i) y F(t ) los vectores de ACELERACIÓN, VELOCIDAD, 
DESPLAZAMIENTO Y FUERZA 







x i 


>1 

x (0 = 


; x (t) = 


; x (t) = 


1! 



*2 . 


*2. 


X 2 . 


F 2. 


Un ejemplo de sistemas con tres grados de libertad es el pórtico simple 
indicado en la Figura 2.3.. Considerando que los dinteles son infinitamente rígidos 
en relación a los pilares, se introduce la hipótesis de que las uniones entre ellos no 
tienen posibilidad de rotación. Por otra parte suponiendo que toda la masa está 
concentrada en los dinteles se transforma un problema con infinitos grados de libertad 
en un sistema con tres grados de libertad, que corresponden a los desplazamientos 
horizontales de las masas. 

De forma análoga a lo realizado en el ejemplo anterior, pueden escribirse las 
ecuaciones del movimiento 

m i y\ = F i + *2 (y 2 - y t ) - K y x 

m 2 y 2 -F 2+ k,(y,- y 2 ) - k, (y 2 - y¿ 
m 3 y 3 = F 3~ k 3 3 - yj 


n.5 



o bien 

m l y\ + (k t + y x - k 2 y 2 = F { 

m i Si + ( k 2 + y% ~ K y\ -^ 3 = F 3 

m 3 y 3 + k 3 y 3 - k 3 y 1 = F 3 

y„(t> 



Fig.2.3. 


De nuevo, se obtiene un sistema de ecuaciones diferenciales acopladas, que 
en forma matricial se escribe 
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A la vista de lo anterior se desprende que el número total de g.d.l. de un 
sistema dinámico será igual al producto del número de masas por el número de tipos 
diferentes de movimientos asociados a cada masa; y que por tanto, la ecuación (2.2-3) 
representa las ecuaciones de un sistema con N g.d.l.. En este caso las matrices son 
de dimensiones (N x N) y las de los vectores (N x 1). Debe observarse que las 
matrices son simétricas y su construcción es automática si se aplican las ideas del 
método directo de la rigidez. 


2.3.- SISTEMAS CON DOS GRADOS DE LIBERTAD 


Al igual que se hizo en el estudio de sistemas con un g.d.l., se comenzará 
resolviendo el problema correspondiente a un sistema sin amortiguamiento y en 
ausencia de solicitaciones exteriores. 

Para ello se considerará un sistema con dos g.d.l.. En las Figuras 2.4. y 2.5., 
se representan dos modelos que permitirían estudiar el movimiento vertical y 
horizontal, respectivamente, de dos sistemas con dos g.d.l.. 



Fig. 2.4. 









k 1 <*1 2* 



Fig. 2.5. 


Si se hace k = O, en el modelo de la Figura 2.4.; las ecuaciones del 
movimiento, en ambos casos, serán 


m i y\ = - k i <yi - 

y 2 = - k 2 y 2 + K (y t - y 2 ) 

o bien 

w i y\ + k i toi - y¿ = 0 

m 2 h + *2 - *1 (y t - yj = 0 

y reordenando 

m i y\ + k i yi - k i y 2 = 0 

m 2 % - *i y i + ( k i + y i = 0 


(2.3-1) 
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(2.3-2) 





Si se buscan soluciones del tipo 


y, = 4>i 


y t " 4> 2 


(2.3-3) 


entonces 

y] = -o 2 4>, e iwt 
y 2 = -u 2 <J> 2 e iwt 


(2.3-4) 


la sustitución de (2.3-3) y (2.3-4) en (2.3-2) proporciona al dividir por e !wt 


-u 2 (j)j nty + ({íj ¿i - <j) 2 = 0 

-<n 2 <}) 2 /n 2 - <j)j + (fcj + ^) <(> 2 = 0 


(2.3-5) 


y reagrupando términos 

- <*> 2 m l ] - <j> 2 ^ = 0 
<J> 2 [ (fc, + k¿ - w 2 m 2 ] - 4>j k x = 0 

las ecuaciones (2.3-6) también pueden escribirse en la forma 
4>i _ (fc, + k 2 ) - a) 2 m 2 

$2 - W 2 ffíj 


(2.3-6) 


(2.3-7) 


lo que implica 

k\ = [ ^ - w 2 /Mj ] [ (fcj + ^) - w 2 m 2 ] (2.3-8) 

la expresión anterior representa una ecuación de segundo grado en to 2 

m i m 2 0)4 “ C m i (^i + & 2 ) + w 2 K í w2 + K = 0 (2.3-9) 

cuyas soluciones son 


w 2 = [ gi + + w 2 ^ ] ± yfr m i ( ^ ^ ) + k, m 2 ] 2 - 4^ m x m 2 

2 m 1 m 2 


(2.3-10) 
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es decir 


id 


[ m l (k x + k 2 ) + m 2 k { ] ± ^[- m l ( k t + k 2 ) + k r m 2 ] 2 + 4 k\ tn l m 2 


2 /Mj m 2 


(2.3-11) 


como puede observarse la cantidad subradical es esencialmente positiva por lo que las 
dos raíces de co 2 , (oí, 2 < o^ 1 ), son reales, denominándose frecuencias propias del 
sistema. La mas baja, Wj, es la frecuencia fundamental y la siguiente, o^, primer 
armónico. 


Si por simplificar el ejercicio se toma 


/Mj = m 2 = m 
k x = k 
&2 = 2 k 


la ecuación (2.3-11) se reduce a 


0) 


2 



\j&k 2 m 2 


= - (2 ± y/2) 

m 


(2.3-12) 


<*>! = — (2-^/2) 
m 

<4 = - ( 2 + ^ ) 

m 


(2.3-13) 


Por otra parte, teniendo en cuenta las dos raíces de to 2 , la ecuación (2.3-7) 
puede escribirse como 


o bien 


4>Í 


<|>2 k x - 4 m { 


2 

(k { + & 2 ) - (x)¡ m z 


(2.3-14) 


¿ W i m i , i . . 2 m i , i 

<j, 2 = <J>1 = ( 1 " «i — ) ^ 

K l 


(2.3-15) 
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como se ha tomado ni, = m y k x = le, (2.3-15) se transforma en 


44 = ( 1 - «? v ) d>i (2.3-16) 

k 

si la componente <¡> x l se hace igual a la unidad, entonces 

<4 = c i - j «í ) = -i + 

(2.3-17) 

Ü = ( 1 - ^ «2 ) = "I - 

k 


A los vectores 


1 


i 


; <l> 2 = 


. -><72 . 


-1-V2 


(2. 3- 17b) 


se les denomina MODOS DE VIBRACIÓN. Puede observarse, Figura 2.6., que en 
el modo fundamental, el movimiento de las dos masas se produce en el mismo 
sentido; mientras que en el segundo modo de vibración, correspondiente al primer 
armónico, el movimiento está en oposición de fase. 


Modo Fundamental de Vibración 





2 ‘ Modo de Vibración 



Fig. 2.6. 
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El planteamiento anterior también podría haberse realizado en forma matricial; 
en efecto, en esta forma la ecuación (2.3-2) puede escribirse como 


m y\t) 

+ k y(t) = 0 


(2.3-18) 

o lo que es igual 

+ 

o 

s 

L— 

k x -ky 

' y» o ' 

(2.3-19) 

S ?s 

1 

5 

O 

-ky ky+k 2 

1 

o 



y buscando soluciones del tipo 

y(f) = 4> e lat (2.3-20) 


se obtendría para el vector de aceleraciones 

y ( i ) = -ü> 2 <|) e iwt (2.3-21) 

la sustitución de (2.3-20) y (2.3-21) en (2.3-1B) proporciona 

-tú 2 m 4> + k 4> = 0 (2.3-22) 

o bien 

(k - Q 2 m) 4> = 0 (2.3-23) 



puesto que se trata de un sistema de ecuaciones homogéneo, la condición de 
existencia de una solución distinta de la trivial se obtendrá al igualar a cero el 
determinante de la matriz de coeficientes del sistema 
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fcj - O) 2 771, -fc, 

-k t ( ^+¿2 ) - w 2 m 2 


= O 


(2.3-26) 


es decir 


[ - <o 2 tu, ] [ (&, + A^) - (o 2 wi 2 ] - kf = O (2.3-27) 


que representa la misma condición expresada por la ecuación (2.3-8). 


Una importante propiedad se obtiene al hacer los productos escalares de los 
vectores (2.3-17b), ponderados por la matriz de masa y rigidez. En efecto, si 






(2.3-28) 


estos productos escalares ponderados por la matriz de masa proporcionan 



4>“ 





<J> ir m cj) 1 

4> U 771 <j) 2 

M = <J) Í m <J> = 

4> 2< 

m 

4» 1 

<j> 2 ' 

= 

<J) 2< 771 (J) 1 

- - 

4> 2< 771 <t> 2 


(2.3-29) 


de forma análoga para la matriz de rigidez, se obtendría 







<j> u k 4) 1 

4) 1 * k (J) 2 

1 

1 

1 

1 

k 

4> l 

4> 2 





<t> 2 ‘ 



<j) 2< k (j) 1 

<J> 2 ‘ k 4> 2 


(2.3-30) 


Para este caso particular, teniendo en cuenta la simplificación realizada en 
(2.3-12), las matrices de rigidez y de masa tomarían la forma 


II 

1 

-1 

J 771 = 771 

1 

0 


-1 

3 

- 

0 

1 


(2.3-31) 
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como 


<t>‘ = 

1 

; 4> 2 = 

1 


-l+l/5 


-1-1/2 


entonces 

4» 1 ' = [1 ; -l +v /2] 

<t> u = [1 ; -1-/2] 

o lo que es igual 


4» = 

1 1 

1 

; 4>' = 


+ 

i-H 

1 

-l-l/5 



si se realizan los productos 

se obtendría 

<{> u m 4) 1 = [ 1 

<j) lf m <J> 2 = [ 1 

4> 2Í m (j) 1 = [ -1 
<|> 2Í m (j) 2 = [ -1 


<J) ir m «J) 1 


(j) 1( m <J) 2 


<j) 2í m (J) 1 


<j) 2t m <J> 2 




-1+1/2] 


m 


m 


1 0 
0 1 

1 0 
0 1 


1 

-l+l/2 

1 

-l-y/2, 


-W5] 


-i +i/5] 


/M 


/n 


1 0 
0 1 

1 P 

0 1 


1 

-1+V5 

1 

-l-l/5 


(2.3-32) 

(2.3-33) 

-l+l/5 

-l-i/2 

(2.3-34) 

(2.3-35) 

= 2/n (2-y^) 

= 0 

= 0 

= 2m (2+v/5) 

(2.3-36) 
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y para la matriz de rigidez 


-l+v/2 ] k 

1 -1 

1 

-1 3 

-l + v/2 . 

-l+v/2 ¡ k 

1 -1 

1 


-1 3 


-l+v/2 ] k 

1 -1 

1 


-1 3 


-1+fi ] k 

1 -1 

1 

-l-v/2 

- 1 3 


= 4 k (3 -2^/2) 


= O 


= O 


= 4 k (3 +2^/2) 


(2.3-37) 

Los resultados obtenidos en (2.3-36) y (2.3-37) indican que, en estos 
particulares productos escalares los modos de vibración son ortogonales respecto de 
la matriz de masa y rigidez. Esta propiedad es muy útil, pues permite un cambio de 

coordenadas a la base 4> que desacopla el sistema de ecuaciones diferenciales que 
caracteriza el comportamiento del sistema. En efecto, si se hace 

y = Z (2.3-38) 


es decir 


o bien 


y 2 



v 

O, 2 


«m 

*-* 


4*2* 

4>2 2 


. ^ . 


yi 


4> l 

4> 2 1 

' «i 

. y 2 . 



- J 

. *2 . 


(2.3-39) 


(2.3-40) 
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(2.3-41) 


que también puede escribirse como 


y = t }) 1 + ( j ) 2 \ 2 


por lo que el vector de aceleraciones, quedaría 

y = 4» 1 I x + l 2 


sustituyendo (2.3-41) y (2.3-42) en (2.3-18), se obtiene 

m 4» 1 + m <t> 2 Í 2 + k 4» 1 ^ + k ( j) 2 í 2 = 0 (2.3-42) 


Premultiplicando (2.3-43) por 4> lf y haciendo uso de (2.3-36) y (2.3-37), se 

obtiene 


(<t> lf m (J) 1 ) fj + (<j> lf * 4» 1 ) ^ = 0 (2.3-44) 


análogamente, premultiplicando (2.3-43) por <|> 2f 

(4> 2 ' m 4> 2 ) l 2 + (<J> 2 ' * <|) 2 ) = 0 


(2.3-45) 


Por otra parte, como, véase (2.3-22), 

-o) 2 m 4> + k <j> = 0 


la premultiplicación de la expresión anterior por 4> ir y <j> 2t , proporciona 


(<t> lf k <j>*) = o), 2 (4> 11 m (j) 1 ) 


(<J> 2 ' k 4> 2 ') = w 2 2 ((J) 2Í m <f> 2 ) 


(2.3-46) 
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en consecuencia, las expresiones (2.3-44) y (2.3-45), pueden escribirse en la forma 

U + “i 2 = ° 


(2.3-47) 


U + W 2 2 *2 = 0 


que son las ecuaciones tipo de sistemas con un g.d.l., cuya solución era 

Ijj = A 1 eos (o>j t) + 2?j sen (o^ t) 

Z 2 = A 2 eos (<i> 2 t) + B 2 sen (w 2 t) 


(2.3-48) 


sustituyendo (2.3-48) en (2.3-41), resulta 


?! 



1 4>i 2 

1 - 


A { eos (Wj t) + fij sen ((o x t) 

. ^ . 


•O-* ! 

1 <j> 2 2 


A 2 eos (o> 2 í) + B 2 sen (o > 2 r) 


<j>{ [ A { eos Wj t + fij sen <o 2 t ] + <J)j [ A 2 eos « 2 t + J3 2 sen o> 2 í J 


[4>2 [ \ eos w x t + B x sen co 2 t ] + 4> 2 [ A 2 eos <ú 2 t + B 2 sen « 2 t ]J 

(2.3-49) 

las constantes A,, B¡, A 2 , B 2 , se determinarán al imponer las condiciones iniciales en 
movimientos y velocidad. 

Como 


y = <J >5 ; y=<M 


(2.3-50) 


entonces, si e yo representan las condiciones iniciales en movimientos y 
velocidad, se tendrá 

y 0 = 4> ; y 0 = <f> Z 0 (2.3-51) 
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premultiplicando (2.3-51) por 4> í m 

(4>' ni) y a = (4>' m <f>) 

(<J)' ni) y a = (4> f m 4>) t o 

por lo que es inmediato escribir 

Z 0 = AT 1 (4)' ni) y a 
= AT 1 (4>* ni) y o 

siendo M 1 i a inversa de la matriz definida en (2.3-29). 
En este ejemplo 


(2.3-52) 


(2.3-53) 


M = 2m 

2-/2 

0 

; A/ 1 = — 

2+/2 

0 


0 

2+/2 

4 m 

0 

2-/2 


obteniéndose, para el producto M 1 ( 4> f m ) 


M 1 4>* m = 


2+/2 0 

0 2-/2 J 

2+i/2 s/2 

2-/2 ~/2 J 

y si, por ejemplo, las condiciones iniciales son 


y 0 = 


1 v/2-1 

1 1 -V5-1 J 


m 


1 0 
0 1 


J. 

4 


• Vl 0= 0) 


5 


y 1 0=0) 


0 

(f=0) 


0 

; y 0 =■ 

y, 

2 0=0) 


0 
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teniendo en cuenta (2.3-53), se obtiene 



2+v/2 

& 

' 6 

_ _6 

2 +v /2 

_ 4 

2-v^ 

-0. 

0 

4 

2-72 


í 0 


por otra parte, como 


resulta 


para t=0 


Si = A t eos (q, t) + B l sen (co 1 t) 
i 2 = A 2 eos (o 2 t) + B 2 sen (co 2 t) 

= “i ( ~ A í sen (cOj t) + B t eos (g)j t) ) 

U = W 2 ( ~ Á 2 Sen («2 *) + B 2 COS («2 l ) > 


Si = A i 
1 (<= 0 ) 1 




(t=0) 


= A „ 


es decir 


= ojj B t 

í 2 = B 2 


L = 


; = 


w i B i 
C0 2 


la identificación de la expresión anterior con (2.3-54), permite obtener 

A i = 6 ; B, = 0 

4 

. í . 0 

4 ^ 


(2.3-54) 


(2.3-55) 
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en consecuencia, de acuerdo con (2.3-48) 


? = 5 hfe. eos ¿ 

1 4 1 


U = 8 


2-V2 


COS 0),í 


Por lo que, finalmente 
5 


= — [ (2+^/2) COS Qjí + (2-y/2) COS ü> 2 f ] 


>-2 


_ 6 y/2 


[ COS Oó.f - COS (0 2 r ] 


2.4.- FRECUENCIAS PROPIAS Y MODOS DE VIBRACIÓN 

La solución del sistema de ecuaciones diferenciales acopladas (2. 2. -3), que 
caracteriza el comportamiento dinámico de un sistema, se basa en las propiedades del 

problema de vibraciones libres, correspondiente al caso c = 0 , F = 0 , por lo que 

la ecuación resultante es la (2.3-18) 

m y(t) + k y(t) = 0 (2.4-1) 

donde, las dimensiones de las matrices, de rigidez y de masa, son (NxN), y las de 
los vectores de movimientos y aceleración (Nxl). 

Con soluciones del tipo 

y(t) = <{) (2.4-2) 


donde 4* no depende del tiempo; se está representando un movimiento armónico de 

pulsación w, tal que todos los puntos pasan simultáneamente por sus valores máximos 
y mínimos. 
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El vector de aceleraciones, se obtendría al derivar dos veces, con respecto al 
tiempo, (2.4-2) 


y (t) = -o) 2 <J) e iM 

la sustitución de (2.4-3) y (2.4-2) en (2.4-1) proporciona 

-o» 2 m <j> + k <J> = 0 


(2.4-3) 


(2.4-4) 


o bien 


(k - o) 2 m) <{> = 0 (2.4-5) 

Por otra parte la ecuación (2.4-5), podría escribirse también en la forma 

k 4> = o) 2 m <J> (2.4-6) 

que representa un problema de autovalores generalizado. En efecto, si se 
premultiplica (2.4-6) por la inversa de la matriz de masa, y se designa por / a la 

matriz identidad de orden N, se tendría 

( m~ l k ) <t> = w 2 / 4> (2.4-7) 


o lo que es igual 


[(m 1 ^)-Q 2 /]4)=0 


(2.4-8) 


(2.4-8) corresponde a un problema de autovalores típico, donde w 2 y <|) son los 
autovalores y autovectores asociados a la matriz A = m 1 k . 
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Los autovalores se obtienen al imponer la condición 


k - m 2 m 


= 0 


(2.4-9) 


que es una ecuación de grado N en o> 2 . Las N soluciones serán reales, puesto 
que k y m lo son. Estas N raíces 

Wj < w 2 < ü> 3 < ... < <¿ N (2.4-10) 

se ordenan de menor a mayor. La mas baja «j, se denomina frecuencia fundamental, 
y a las restantes primero, segundo, etc. armónicos. 

Al sustituir cada valor de oj¡ en (2.4-5) es posible determinar el autovector 

correspondiente <t)‘ , tal que 

k <}>* = g> 2 m 4>‘ (2.4-11) 


obteniéndose de esta forma los modos de vibración. 

Debido al carácter homogéneo del sistema de ecuaciones (2.4-5), los modos 
están indeterminados por lo que se puede hacer la primera componente de cada modo, 

<J>‘ , igual a la unidad, es decir 



(2.4-12) 


aunque en otros caso se opta por normalizar el vector 

£ i V f - i (2 - 4 - ]3) 

¡t= i 

o dividir las componentes por el valor máximo de las mismas. También, con objeto 
de simplificar las formulas, algunos autores patrocinan una normalización del tipo 

<}>* ‘ m <}>* = 1 (2.4-14) 
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EJEMPLO 2.1 


Obtener las frecuencias propias y modos de vibración del pórtico indicado 
en la Figura E.2.1-1. Los pilares se supondrán inextensibles y sin masa ; los dinteles 
se considerarán infinitamente rígidos en relación a los pilares y con masas m, 1.5m 
y 2m. 


Datos : El — 8 x 10 ( N x m 2 ; m = 5000 kg ; L = 4 m 



Fig. E.2.1-1 


Solución : 
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* La matriz de rigidez de la estructura, para los grados de libertad indicados se 
determinará teniendo en cuenta que la columna i son las fuerzas que aparecen en los 
g.d.l. j al dar un desplazamiento unidad en el g.d.l. i, manteniendo impedidos el resto 
de los grados de libertad. De este modo 

- I a Columna : kj, , j = 1, 2 y 3 


* 1 tK * 1 tK k lt 



- 2 a Columna : k j2 , j = 1, 2 y 3 

tK * -X- 7^ 1 -)f *12 



- 3 a Columna : kp , j = 1, 2 y 3 

* 1 * tK 1 # k 13 
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NOTA : Obsérvese que los coeficientes han sido obtenidos teniendo en cuenta que 


6EI 6 6EI 6 



por tanto 


24 £7 
L 3 


1 

-1 

0 


-1 

3 

-2 


0 

-2 

5 


siendo la matriz de masa 


m = m 


1 

0 

0 


0 o 

1.5 0 

0 2 


Las frecuencias propias se obtendrán al imponer la condición (2.4-9) 

I k - (o 2 m I = 0 


que para este ejemplo se transforma en 


24£7 

£ 3 


1 

-1 

0 


-1 

3 

-2 


0 

-2 

5 


co 2 m 


1 

0 

0 


0 

1.5 

0 



(E2.1-1) 


o bien al dividir por 24EI/L 3 

1 -1 0 

-1 3 -2 

0-2 5 


w 2 m 


Ü 


24EI 


1 

0 

0 


0 

1.5 

0 



(E2.1-2) 
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y si se designa por 


a = ü > 2 m 

24EI 

la condición (E2. 1-2) puede expresarse como 

1-a -1 0 

-1 3-1.5a -2 

0 -2 5 -2a 


= 0 


que proporciona una ecuación cubica en a 

a 3 - 5.5a 2 + 7.5a -2=0 

cuyas soluciones son 

a l = 0.351 
a 2 = 1.608 
a 3 = 3.541 

entonces, según (E2.1-3) 


2 a 24EI 
ar = 


m l 3 


por lo que las frecuencias propias serán 


Qj = 2.9024 


(ú 2 = 6.2122 


El 

mÜ 


El 

mi 3 


rdfseg ; T x 


rdfseg ; 7, 


o) 3 = 9.2187 


N 


El 

mÜ 


rdjseg ; 7, 


= 0.433 seg 

= 0.2023 seg 

= 0.1363 seg 


(E2.1-3) 


(E2.1-4) 


(E2.1-5) 


(E2.1-6) 


(E2.1-7) 
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En lo que a los modos de vibración se refiere, éstos se obtendrán imponiendo 
la condición 


(k - w 2 m) <f > ¡ = 0 ; j = 1, 2, 3 

o lo que es igual 


1-a^ -1 0 


44 


0 

-1 3-1. 5a. -2 


44 

= 

0 

0 -2 5-2a. 


tj 


0 

i 


<P 3 




y haciendo la primera componente de cada modo igual a la unidad, es decir 


= i - <i>l = 4>í = <t>í = i 

se obtiene 


1-CCy -1 0 


1 


0 

-1 3-1.5a^ -2 


44 

= 

0 

0 -2 5-2ctj 


44 


0 . 


de la primera de las ecuaciones 

1 - = 4 4 

y de la tercera 

2 <|4 = (5-2a.) 44 

por tanto 

44 = I-a, 
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en consecuencia, se tendrá 


* Primer Modo de Vibración 
T, = 0.433 seg 

4>J = 1 

<J>Í = 0.649 
<t>3 = 0.302 


* Segundo Modo de Vibración 
T 2 = 0.2023 seg 

<D? = l 

<{>2 = -0.608 

4 >¡ = -0.682 


* Tercer Modo de Vibración : T 3 = 0.1363 seg 

4>; = i 
<j>2 = -2.541 
<t>3 = 2.441 


En la Figura E.2. 1-2 se representan los tres modos de vibración, obtenidos al 
hacer la primera componente de cada modo igual a la unidad. Obsérvese que al estar 
los modos indeterminados cualquiera de las cuatro posibles configuraciones 
representadas en la Figura E.2.1-3 son también aceptables. 





Pri«er Modo 


Segundo Modo 


Tercer Modo 


Fig. E.2.1-2 
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en consecuencia, los Modos de Vibración son 


* Primer Modo 
Tj = 0.433 seg 

♦ 1=1 
<|>‘ = 0.649 
<j>‘ = 0.302 


* Segundo Modo 
T 2 =0.2023 seg 

1 

^ = -0.608 
^ =-0.682 


* Tercer Modo 
T 3 =0.1363 seg 

4>, 3 =1 

♦* = - 2.541 
4>3 = 2.441 


En la figura E.2.1-2, se representan los tres modos de vibración, obtenidos al hacer la 
primera componente de cada modo igual a la unidad. Obsérvese que, al estar los modos 
indeterminados, cualquiera de las cuatro posibles configuraciones representadas, para cada 
modo, en la Figura E.2. 1-3 son también aceptables. 

i < i 





Fig. E.2.1-2 



Fig. E.2.1-3 



4 3 


2 1 



Prleer Modo 


4 3 2 1 



Segundo Modo 


4 3 2 1 



Tercer Modo 


Fig. E.2.1-3 


2.5.- PROPIEDADES DE LOS MODOS DE VIBRACIÓN 


Como cada uno de los modos de vibración ha de cumplir la condición (2.4-6), 
es posible escribir para el modo j 

k = w? m <J) y (2.5-1) 


Por otra parte como k y m son matrices simétricas se verificará 

m t = m ; k‘ = k (2.5-2) 

siendo entonces la traspuesta de (2.5-1) 

($>*)* k = Wf 2 (<j> ¿y m (2.5-3) 
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(2.5-4) 


Para una frecuencia w k w¡ 

k 4>* = u 2 m 4>* 

y post multiplicando (2.5-3) por <J>* 

(4>Ó r k <}>* = w. 2 (4)^ m 4>* (2.5-5) 

sustituyendo en la expresión anterior, k <j>* , por su valor dado en (2.5-4), se obtiene 

<¿ k 2 (4)0' m 4)* = o* (<f) Ó* m 4>* (2.5-6) 

o lo que es igual 

((ú k 2 - <ú?) [ (4>Ó' m <b k ] =0 (2.5-7) 

y como w k coj resulta 

(4>Ó t w4>* = 0 (2.5-8) 

Análogamente, la sustitución de, m 4>* , por su valor dado en (2.5-4), en 
(2.5-5), permite escribir 

<ú k 2 (4)^ k <$> k = tí 2 (fr*)' k <b k (2.5-9) 

es decir 

(<*>/ ~ <•>/) [ (4> / )' k 4>* ] =0 (2.5-10) 
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y como Wk ^ Wj se obtiene 


(2.5-11) 


(<t> ^ k 4>* = 0 


Si las matrices de masa y rigidez fuesen la matriz identidad de orden N, 
entonces 

m = I ; k = / 


I 

(4>^ f . <j>* = o 

f 

<t>' -L <j>* 


en consecuencia (2.5-8) y (2.5-11), en este caso particular, expresan la condición 
clásica de ortogonalidad entre dos vectores. Siendo, por tanto, natural considerar a 
las expresiones 


(<J>Ó f tn <i>* 
(4»^ k (J>* 


como productos escalares generalizados, y afirmar, tal y como se apuntaba en el 
apartado 2.3, que en este producto escalar los modos de vibración son ortogonales 
respecto de la matriz de masa y rigidez. 

En el caso en el que co k = Wjj es decir, cuando k = j, los productos 

(<!>*)' m <j>* 

(4>V i fe 


no se anulan; siendo habitual designar como masa modal al producto 

M k = (<!>*)' m <f>* (2.5-12) 


y normalizar los modos 4>* de manera que se cumplan la condición M k = 1, con lo 
que se simplifica algo la notación. 
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En resumen, las expresiones (2.5-8), (2.5-11) y (2.5-5) implican 
k * j 

(4>^ f m 4>* =0 

(4) Ó r k 4>* = 0 (2.5-13) 

k = i 

(4»*)' k <j)* = g>¿ 2 M k 


Por último, y como consecuencia de lo anterior, cabe indicar que los 
productos 


4»* m <j) ; 4>* k 4> 


son matrices diagonales que son designadas por 

M = & m <J> ; í: = 4>‘ k<b 


(2.5-14) 


siendo 


4> = 


4» 1 I ! ... I <J> W 


4»! 1 4 > i 2 - 4 >/ 

V $2 ••• ^ 

••««••4 444 4444444 

V V ... 4>/ 


(2.5-15) 


la denominada matriz modal, cuyas columnas son los modos de vibración. 
Así pues 

0 ; Vi * k 


M ik = (4>y m 4>* = 


M. ; V i = k 


(2.5-16) 


K ik = (4>‘)' k <j> 


k _ 


0 ; V i * k 
g> 2 M.;Vi = k 


(2.5-17) 
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2.6.- EL MÉTODO DE LOS MODOS NORMALES 


Aprovechando las propiedades de ortogonalidad de los modos y puesto que 
éstos son soluciones del sistema homogéneo surge la idea de uti liz ar una combinación 
de los mismos para calcular la solución general. En otras palabras, se propone 
calcular ésta mediante la combinación lineal 

y = <J> 5 = (j) 1 + 4> 2 $ 2 + ... + 4> n (2.6-1) 


lo que indica claramente un cambio de coordenadas a la base modal. En estas 
condiciones la ecuación (2.4-1) se transforma en 

/M<t>í;+fc<J>£=0 (2.6-2) 


Premultiplicando (2.6-2) por <J) f , se obtiene 

(<!>' m 4») % + (<{>' k <t>) $ =0 (2.6-3) 


teniendo en cuenta las propiedades de ortogonalidad de los modos y (2.5-14), la 
expresión anterior toma la forma 

M l + K £ = 0 (2.6-4) 


al ser las matrices M y K diagonales, véase (2.5-16) y (2.5-17), el sistema (2.6-3) 
representa N ecuaciones diferenciales correspondientes a un sistema con un g.d.l.. 


o 

o 

sf 




GJj 2 Mj 

o 

sT 

o 


u 


0 




+ 


0 0 ... m n 


Ik. 


0 


0 ... 0 


V 


0 

o 

sf 

3 


Í2 


0 





• *« 

0 ... o>„ 2 M n 


% N . 


0 


(2.6-5) 
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Es decir, la resolución del sistema de N ecuaciones diferenciales acopladas, 
se ha transformado en el análisis de N ecuaciones diferenciales del tipo 


%k + l k = 0 


(2.6-6) 


que son las ecuaciones tipo de sistemas con un g.d.l., cuya solución es 


Z k = A k eos (tó* t) + B k sen t) 


(2.6-7) 


las constantes A k y Bj., se determinarán al imponer las condiciones iniciales en 
movimientos y velocidad. 

Tal y como se comentó en el apartado 2.3, si y o e y o representan las 
condiciones iniciales en movimientos y velocidad, se tendrá 

y a = 5 , ; y 0 = 4 > i 0 ( 2 . 6 - 8 ) 


premultiplicando (2.6-7) por ¥ m 

(<J>' m) y 0 = (<{>' m <{>) 
(<!>' m) y o = (4> f m <J>) 


(2.6-9) 


por lo que es inmediato escribir 


l 0 = M 1 (<}>' m) y a 

K = M ~ l (¥ y o 


( 2 . 6 - 10 ) 
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siendo M 1 la inversa de la matriz definida en (2.5-16); es decir 


L = 


M, 


0 


0 


Ai, 


0 0 


0 

0 


M, 


N 


4 >' « y„ 


( 2 . 6 - 11 ) 


= 


— 0 
M, 


0 — 
Ai, 


0 0 .. 


0 

o 


Ai 


N 


4>' m y 0 


( 2 . 6 - 12 ) 


conociendo los valores de % a y Í a , es inmediato obtener las constantes A k y B k , 


ya que, para t=0 


( tk X> = 


(2.6-13) 


por lo que finalmente 


= ( \ k )„ eos (<¿ k t) + - — — sen (<ó k t) 


(2.6-14) 


(ú í 


2.7.- VIBRACIONES FORZADAS 


En este apartado se aplica el método de los Modos Normales al sistema 

my' + ky = F (2.7-1) 
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donde F(t) es el vector de fuerzas en los grados de libertad al discretizar la 
estructura. Utilizando el cambio de base 


y = <i> 5 ; y = <J> l 


(2.7-1) se transforma en 


m 4> S + £<t>£ = F 


premultiplicando (2.7-3) por 4)' , se obtiene 


($* m 4>) l + (<J)' k <{>) £ = 4>* F 


y como 


(4>') í m 4>* = 


0 ; V i * k 
M t ;V i = k 


(4> f )' k 4)* = 


0 ; V i * k 
w 2 ; V i = k 


la ecuación (2.7-4) puede escribirse 


li 


+ <i). 


Ci 


(4>*y F 


ecuación análoga a la obtenida en apartados previos, salvo que 
miembro aparece la proyección del vector de cargas sobre el modo 
Para poner de manifiesto la simetría del proceso de proyección, en 
escribirse el segundo miembro de (2.7-5) en la forma 

(4>y F (4>y F 


(4>y m <J>' 


(2.7-2) 


(2.7-3) 


(2.7-4) 


(2.7-5) 


en el segundo 
de vibración i. 
ocasiones suele 


(2.7-6) 
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La solución de (2.7-5) se calcula en la forma habitual para sistemas con un 
g.d.l., es decir 


= A. eos (g>. i) + B. sen (w ¿ t) + (2.7-7) 

donde £ p ¡ es la solución particular de la ecuación (2.7-5), que puede obtenerse por 
cualquiera de los procedimientos indicados en el Capitulo I (vg.: la integral de 
Duhamel). 

Las constantes A¡ y B¡, se determinarán a partir de las condiciones iniciales 
(£y) 0 y (Íj) 0 , definidas en (2.6-13); en efecto para t=0 (2.7-7) proporciona 

((, )„ - (, (<=0) . A, * (Íí )„ (2-7-8) 

y como 

í¡ = «. [ -A t sen (w j t) + B. eos (w ¿ t) ] + ^ (2.7-9) 

se obtiene 

(*, \ = z¡ < í= o) = <*>, b. + a p i \ < 2 - 7 - 10 ) 


2.8.- RESPUESTA DE SISTEMAS AMORTIGUADOS 


En este caso la ecuación que caracteriza el comportamiento del sistema es 

my+cy+ky=F (2.8-1) 

y realizando el cambio de base 

y = 4>£ ; y = < \> Í ; y = <j> 1; (2.8-2) 

(2.8-1) se transforma en 

m$l+c$Í+k$f,=F (2.8-3) 
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premultiplicando (2.7-3) por 4 / , se obtiene 

( 4 >' m 4>) l + (<})' c 4 >) í + ( 4 >‘ k 4 >) $ = 4 >‘ F (2.8-4) 


mientras que (4> r 4>) y (4>* k 4>) son matrices diagonales; no existe, sin 


embargo, ninguna condición sobre el término (4> r c 4>) , por lo que el 

desacoplamiento no es total. Para conseguirlo es preciso añadir alguna hipótesis sobre 
la matriz de amortiguamiento, siendo la mas inmediata la denominada 
amortiguamiento proporcional, es decir 

c = a m + p k (2.8-5) 


donde a y 0 son parámetros que dependen del tipo de estructura por lo que deberán 
fijarse para cada caso particular. Con la hipótesis previa, premultiplicando (2.8-5) 

por 4> ‘ y postmultiplicando por 4 ) , puede escribirse 

4>' c 4> = a (4>' m 4>) + p (4>' k 4>) (2.8-6) 

o lo que es igual 

4> í c4) = aM + p J &: (2.8-7) 


y al ser Ai y K , matrices diagonales (4> 1 c 40 también lo será; pudiendo escribir 
para la ecuación i-esima 


M, l¡ + (a Ai, + fJ cof Ai.) £, + o), 2 Ai, $, = (4^)' F (2.8-8) 


y dividiendo por M¡ 


f 

+ (a + p a>f) + (o, 2 $, = — - 

Ai; 


(2.8-9) 
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y si se hace 


2( i o ¡ = a + 


donde 



c ¡ 


c i 

l crit 


la ecuación (2.8-9) toma la forma 


l + 2^ É, * 


(*')' F 




(2.8-10) 


( 2 . 8 - 11 ) 


( 2 . 8 - 12 ) 


La relación (2.8-10) puede utilizarse para obtener a y /?, cuando se conocen 
los índices de amortiguamiento para dos modos cualesquiera. Así, por ejemplo, si se 
conocen los índices de amortiguamiento para los dos primeros modos, y 

2 (. 0 l>. = a + Po>, 

(2.8-13) 

2í 2 w 2 = a + pa >2 


que representa un sistema de dos ecuaciones y dos incógnitas, ay/3 


a 


2 g > 1 (») 2 


íl W 2 ~ C 2 <*>i 


6), 


W, 


P 


= 2 


C 2 a > 2 

(0 2 - 


(2.8-14) 


obtenidos los valores de « y /3, el índice de amortiguamiento para un modo cualquiera 
i, se obtendrá teniendo en cuenta (2.8-10) 



a + Pwf 
2w i 


(2.8-15) 
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para el caso particular = j* 


a 


2Í 


<*l <*>2 

6 ), + (ó 2 


P 


2C 


1 

w, + < ó 2 


(2.8-16) 


La solución de (2.8-) se reduce, de nuevo, al caso del sistema con un g.d.l. 

eos (<o^ t) + B i sen (w irf t) + (2.8-17) 


donde 

= (i)¿ ^1 - Cf (2.8-18) 

y es una solución particular que debe calcularse para cada grado de libertad. Las 
condiciones iniciales permiten, a su vez determinar los valores de las constantes A- 

y b¡. 


2.9.- MOVIMIENTO E N LA BASE, VECTOR DE ARRASTRE. FACTOR DE 
PARTICIPACIÓN 

Un caso importante de carga es el originado por un movimiento de la base; 
por ejemplo, el producido por un terremoto 



Fig. 2.7. 
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La carga total que actúa en cada grado de libertad puede obtenerse 
generalizando, la situación expuesta en el Capitulo I para sistemas con un g.d.l., al 
caso de sistemas con N g.d.l. cuya unión con el suelo sufre un desplazamiento y s (t). 

El movimiento total 


y\t) = y s (t ) J + y(t) 


(2.9-1) 


se compone de 

* uno de arrastre 

yj® J (2.9-2) 

* otro relativo 


y(t) 


(2.9-3) 


donde J es el denominado VECTOR DE ARRASTRE que indica los grados de 


libertad afectados por el movimiento de solido rígido. 


La Figura 2.8. recoge los vectores de influencia para una estructura simple 
formada por barras y con g.d.l. de traslación 


í 

7^ L 3 -tK 




1 


0 

J = 

1 

J = 

0 


0 


1 


MOVIMIENTO MOVIMIENTO 
HORIZONTAL VERTICAL 
EN LA BASE EN LA BASE 


Fig. 2.8. 
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Con estas consideraciones la ecuación del movimiento es 


m y T + c y + ky = F (2.9-4) 

y teniendo en cuenta ( 2 . 9 - 1 ) 

m y + cy+ky = -mJ y s (t) ( 2 . 9 . 5 ) 

la ecuación anterior es semejante a la ( 2 . 8 - 1 ) si se hace 

F = -mJ y s (t ) (2.9-6) 


surgiendo la fuerza, en este caso, por efecto de la inercia provocada ante el 
movimiento en la base. 


Utilizando la transformación 

? = 4 > 1; ; y = <M ; y = 


la ecuación (2.9-6) puede escribirse en la forma 


l + 2 < |W| ♦ W| 2 5| 


(<)>')' m J 
m 4 »' y * (t) 


(2.9-7) 


es decir, el estudio del sistema con N g.d.l. se reduce de, esta forma, al de N 
sistemas con un g.d.l. y cargas 


donde 


m = - r ¡ m 


(2.9-8) 


(<{)')' m J 
(<J)'y m (}>' 


(2.9-9) 


es el denominado FACTOR DE PARTICIPACIÓN ya que puede utilizarse como 
criterio para analizar la importancia relativa de los modos de vibración en la 
superposición final. 
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El significado de este factor, se obtiene al desarrollar el vector de influencia 
J en función de los modos. En efecto; supóngase que J puede escribirse en la 

forma 


N 


J - E c* 


(2.9-10) 


sustituyendo (2.9-10) en (2.9-9) 


r. = 


N 


E 

[»“)' M C, *‘l 

*= l 



(2.9-11) 


( 4 >')' m <{)' 


recordando ahora la propiedad de ortogonalidad de los modos respecto de la matriz 
de masa 


(<J>'y m 4>* = 0 ; V k * i 


la expresión (2.9-11) se simplifica 

(4 >y m 4>' 

T. = C. — (2.9-12) 

1 ‘ (4> i ) < m 4> j 


o lo que es igual 

I\ = C, (2-9-13) 


por lo que (2.9-10) puede escribirse 

J - £ T f 4>‘ (2-9-14) 

*=i 

lo que pone de manifiesto que el factor de participación expresa la contribución de 
cada modo al movimiento de arrastre. 



2.10.- ANÁLISIS MODAL ESPECTRAL 

Este procedimiento es uno de los más utilizados en el cálculo sísmico de 
estructuras. Las hipótesis fundamentales que se realizan son las siguientes : 

* El movimiento del suelo , que se aplica a todos los soportes 
simultáneamente, se define en forma de Espectro de Respuesta. 

* El comportamiento de la estructura es lineal , y por tanto se puede 
utilizar la superposición de modos. 

El método consiste en la obtención de las frecuencias propias y modos de 
vibración ( resolución del problema de autovalores) y en el cálculo de la respuesta 
máxima de cada modo de acuerdo con el espectro de respuesta elegido. 

Recordando que la solución de un sistema con un g.d.l. caracterizado por la 
ecuación del tipo 

m y(t ) + c y(t ) + k y(t) = -m y s (t ) (2.10-1) 

o bien 

y(t) + 2Cw y(t) + a 2 y(t) = - y/t) 
podía escribirse en la forma 

t 

y(t ) = — — f e “ T) yfx) sen a> d ( t - x) dx 
w d o 

comparando (2. 10-2) con 

Si + Si + “, 2 €, - - r, m 

y teniendo en cuenta (2.10-3), se obtiene para £¡ 

r. * 

Zft) = f e í<w<íí " T) yfx) sen ( t - x) dx 

w d, o 1 

debido a que los índices de amortiguamiento utilizados en ingeniería estructural son 
muy bajos, se puede suponer que w ¡ = o> d¡ , y por tanto 

- — PSV( a ,, (,) (2.10-6) 

0 ). 


( 2 . 10 - 2 ) 


(2.10-3) 


(2.10-4) 


(2.10-5) 
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y como 


PSA = o. PSV (2.10-7) 

resulta que 

^ = — 2 psa ^ í>Q (2.10-8) 

ahora bien, según la transformación 

y - í; ♦* «, - e y < 2 - 10 - 9 > 

1=1 - i=l “ 

el vector de desplazamientos máximos, para el modo i, teniendo en cuenta (2.10-8), 
será 

/ L = 4>' )max = V «A (O),. , Q (2.10-10) 

" G)¡ 


Por otra parte como 

= K y* )*** = m y 1 ^ (2.10-11) 

es decir las fuerzas equivalentes asociadas al modo i, teniendo en cuenta (2.10-8), 
serán 

= «? * ♦* )mM = « ♦* r ¿ «A (“i . Q (2.10-12) 


La combinación de los distintos modos puede realizarse utilizando diversas 
reglas. La más pesimista, que por serlo en exceso no suele utilizarse, es la que 
obtiene la respuesta como suma de los valores absolutos de las de cada modo 

- E 1 z ¡ 1 (2.10-13) 
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La más utilizada en la práctica es la combinación cuadrática (RMS), en la que 
la respuesta es la raíz cuadrada de la suma de los cuadrados de la respuesta de cada 
modo 


"“total 


E 


(2.10-14) 


Esta combinación está basada en la hipótesis de que la respuesta a cada modo 
de vibración es una variable aleatoria gaussiana independiente, de media nula y 
varianza igual a la respuesta máxima. En general, la hipótesis de independencia será 
fiable cuando los modos de vibración estén suficientemente desacoplados ( las 
frecuencias naturales sean muy distintas). Cuando los modos están acoplados, existen 
reglas mixtas que proponen combinaciones especiales para dichos modos. Dos 
ejemplos de este tipo de combinación son : 

* Regla de 10% (TPM). Si la diferencia entre las frecuencias de dos 
modos es menor del 10% , se suma el valor absoluto de la respuesta de 
ambos modos, y el resultado se combina con los demás mediante la 
regla RMS 


< o.i - Z* = i Z. I + I Z. I (2.10-15) 


Regla CQC (Combinación Cuadrática Completa). Propone un 
coeficiente de acoplamiento entre cada dos modos, que depende del 
cociente entre las frecuencias y del índice de amortiguamiento 


"“total 


1 


EEs-iz.iiz, 

i j 


= 1 

(2.10-16) 

<0. 



A pesar de que el Análisis Modal Espectral está estrictamente limitado al 
estudio del comportamiento lineal, debido a la utilización de la superposición modal, 
se suele aceptar la simulación del comportamiento no lineal, introduciendo un factor 
de ductilidad que minora el espectro de respuesta. 
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Una de las ventajas del método de superposición modal consiste en la 
posibilidad de truncar dicha superposición a un número de modos M<N. Es decir, 
si M es el número de modos que van a utilizarse, se trata de que la ecuación 

y - E *' «. 

¿= 1 - 

pueda ponerse en la forma 

y - í: ♦' e, ♦ E ¥ij (2 - 10 - 17) 

¿=1 - j=M+l - 


de este modo no sería preciso el cálculo de los modos altos. El principal problema 
se presenta a la hora de decidir cuántos modos se deben considerar para obtener una 
buena aproximación. 


2.11- MASA MOVILIZADA. NUMERO DE MODOS A SUPERPONER 


La masa total movilizada por el vector de arrastre será 


M r = J* m J 

1 MOV 


( 2 . 11 - 1 ) 


y como, según (2.9-14), se tiene 


N 


J - E r» *>' 


*= 1 


sustituyendo la expresión anterior en (2.11-1), se obtiene 


N 


N 


M. 


r MOV 


■(Er.irti.ijr,* 1 ) 


*= 1 




( 2 . 11 - 2 ) 


realizando los productos y utilizando, de nuevo, la propiedad de ortogonalidad de los 
modos, resulta que la expresión anterior es igual a 


Mr uov - E rí «’*)' *.♦*-£ «, 


*=1 


Jfc=l 


MOV 


(2.11-3) 
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donde 


= d (**)' »• +* 


'■MOV 


(2.11-4) 


puede considerarse como una masa movilizada en el modo j. Teniendo en cuenta (2.9- 
9), la expresión anterior queda como sigue 


(<j>*)' m J 


M, 


= r L 


uov " 1 rn <|>* 
- r, f (♦'y m j 


(<{)*)' m 4>* = 


(2.11-5) 


Utilizando como criterio para el truncamiento modal, el que los M modos a 
superponer movilicen el 90% de la masa total, se tendría 


T M. z 

k w 


0.90 


1 MOV 


( 2 . 11 - 6 ) 


No obstante, es preciso señalar que, en determinados problemas, es importante 
el estudio de la influencia de Modos Locales de vibración, entendiendo por éstos, los 
modos altos que puedan estar incluidos en aquellos que movilizan menos de un 10% 
de la masa total del sistema, y que, por tanto, siguiendo el criterio de truncamiento 
modal, generalmente aceptado, no serían tenidos en cuenta pudiendo, sin embargo, 
incrementar los esfuerzos de forma notable. 


EJEMPLO 2.2 


Calcular la masa movilizada en cada modo de vibración del pórtico del 
ejemplo 2.1.. 


Solución : 
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Los resultados obtenidos en el ejemplo 2. 1 , proporcionaban para los Modos 
de Vibración los siguientes valores 

Primer Modo Segundo Modo Tercer Modo 



1 


1 


1 

4) 1 = 

0.649 

4> 2 = 

-0.608 

4> 3 = 

-2.541 


0.302 


-0.682 


2.441 


a partir de ellos y teniendo en cuenta que la matriz de masa era 


m = m 


1 

0 

0 


0 0 

1.5 0 

0 2 


es posible determinar la matriz de masa modal, véase (2.5-12), 


M = 4> í m <J) = 



1 

0.649 

0.302 


1 

0 

0 


1 

1 

1 

m 

1 

-0.608 

-0.682 


0 

1.5 

0 


0.649 

-0.608 

-2.541 


1 

-2.541 

2.441 


0 

0 

2 


0.302 

-0.682 

2.441 



1 

0.974 

0.604 


1 

1 

1 

m 

1 

-0.912 

-1.364 


0.649 

-0.608 

-2.541 


1 

-3.812 

4.882 


0.302 

-0.682 

2.441 


1.82 

0 

0 





m 

0 

2.49 

0 






0 

0 







n.49 



En el supuesto de un movimiento horizontal en la base, éste afecta a los tres 
grados de libertad definidos por lo que el vector de arrastre será 



1 

1 

1 


y según (2.11-1), la Masa Total 

M r = J‘ m J = 
t mov 


= rn 


1 


1 


1 



1 

0 

0 


0 

1.5 

0 


0 

0 

2 
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= m (1+1. 5+2) = 4.5 m = 4.5 x 5000 = 22500 kg 


Con los cálculos previamente realizados es posible determinar la masa 
movilizada por cada modo 

* Primer Modo 


(4) 1 )* m J = m [ 1 0.974 0.604 


1 

] 1 
1 


= 2.58 m 


= (<f> l ) r m «(i 1 = 1.82 m 


(4) 1 )' m J 


r i = 


M, 


2.58 m 
1.82 m 


1.42 
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M. = r. («J) 1 )' m J = 1.42 x 2.58 m = 

* MOV 1 _ 

= 1.42 x 2.58 x 5000 = 18318 kg 


1 


M, 


1 MOV 


= 81.4 % M t 

t mov 


* Segundo Modo 


(<j> 2 )‘ m J = m [ 1 -0.912 -1.364 ] 


1 

1 

1 


-1.28 m 


M 2 = C^ 2 )* m <{> 2 = 2.49 m 


(<t> 2 ) r m J 


r 2 = 




-1,28 m 
2.49 m 


= -0.51 


M %mov = r 2 ((j) 2 y m J = (-0.51) X (-1.28) m = 


= 0.51 x 1.28 x 5000 = 3264 kg 




‘■MOV 


14.5 % Af T 

t mov 


* Tercer Modo 


(<j) 3 )' mJ = m[ 1 -3.812 4.882 ] 


1 

1 

1 
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M -i = (<J> 3 )' m 4» 3 = m 



(4> 3 y m J 


AL 


-2.07 m 

22.6 m 


= 0.09 


M z mov = r 3 (<|» 3 )' w / = 0.09 * 2.07 m = 
= 0.09 * 2.07 x 5000 = 931.5 kg 


AL 


* MOV 


= 4.1 % 


' MOV 


2.12.- SISTEMAS CON PARÁMETROS DISTRIRI linos 


En los apartados anteriores se ha analizado un procedimiento cuyo desarrollo 
ha sido posible gracias al hecho de haber considerado que las masas que intervenían 
en el problema estaban concentradas en los grados de libertad de la estructura. El 
problema que se plantea ahora es como determinar una matriz de masas discretas 
cuando la masa está repartida de forma continua a lo largo de la estructura. El 
procedimiento que se utiliza es la reducción del número de grados de libertad, es 
decir la discretización de la estructura. Este proceso se apoya en dos puntos 
fundamentales: 

* El Principio de los Trabajos Virtuales (P.T.V.), como expresión del 
equilibrio entre las energías cinética y potencial del sistema. 

* Una aproximación de la solución como combinación lineal de funciones 
conocidas. 

La idea tuvo su origen en Rayleigh, fue continuada por Ritz y su desarrollo 
actual mas conocido es el Método de los Elementos Finitos. 
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En particular, por ejemplo, para el caso de la barra a flexión 



Fig. 2.9. 


si: 

E es modulo de Young del material 
I el momento de inercia de la sección transversal 
q la carga repartida por unidad de longitud 

ipi (x) una función conocida, perteneciente a una familia { \í-¡ (x) } que cumple 
las condiciones de contorno esenciales. Por ejemplo, sen 2xx/L en 
{ sen ixx /L } 

el P.T.V. se escribe 


*, Q 


dty. L 
— M i + 
dx o 


f 


d 2 \ Jj. 

—22 El 

dx 2 dx 2 


cPy 


dx = J q íjr. dx (2.12-1) 


donde los diferentes términos representan 


t Q 


— M 
dx 


: Trabajo de los cortantes 
: Trabajo de los Flectores 


L f El fy & 

o dx 2 dx 2 


Trabajo de los Flectores en las curvaturas ficticias 


L 

J q ijr dx 


O 


: Trabajo de las cargas exteriores 


( 2 . 12 - 2 ) 
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y al cumplir las funciones \p t (x) las condiciones de contorno esenciales, entonces 
(2.12-1) puede escribirse en la forma 


/ 


d 2 ^¡ ¿py L r 

A E ' * = / q *< * 


dx 2 


(2.12-3) 


En el caso dinámico, la carga repartida q está formada por 

* la carga real actuando sobre la pieza p(x,t) 

* las fuerzas disipativas c(x) y(x,t ) 

* las fuerzas de inercia m(x) y(x,t ) 

q(x,t) = p(x,t) - m(x) y(x,t ) - c(x) y(x,t ) (2.12-4) 


donde se ha supuesto que, m(x) es la masa repartida por unidad de longitud, 
dependiente del espacio pero no del tiempo; y c(x) el amortiguamiento viscoso 
repartido por unidad de longitud. 


Sustituyendo (2.12-4) en (2.12-3), y reordenando se obtiene 

/ m(x) ty¡(x) — dx + [ c(jc) ijr (x) dx + 

o dt 2 J 0 dt 

'f^En^dx-fpMWd* 

o dxr dx 2 i 


(2.12-5) 


si ahora se aproxima la solución real mediante una combinación lineal de funciones 
{ 'Pj (x) }, es decir 


N 


y(x,t) * Y a f!) 
n 


( 2 . 12 - 6 ) 


es inmediato 


N 


y(x,t) » Y á fi> 
/=! 


N 


(2.12-7) 


y(x,t) * Y i|r/x) 
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Puesto que (2. 12-5) expresa el reparto de energía, introducida por las cargas 
p(x,t) en la estructura, entre los términos de inercia, disipativo y elástico; surge la 
idea de determinar los N coeficientes aj (t) imponiendo repetidamente (2. 12-5) para 
N funciones (x). Obteniéndose 

N 

E [ ,n ij á j <0 + c íj á j (*) + k a a j (0 ] = Pi (0 (2.12-8) 

V i = 1, 2, , ... , N 

donde se ha hecho 


L 

m ü = / m(x) t| t/x) dx 
o 

L 

C ij = / 'W c (*) %<*) 

0 



L 

Pi (0 = / p(x>t) 'Jí.W dx 

o 


(2.12-9) 


( 2 . 12 - 10 ) 


por lo que en notación matricial, la ecuación (2.12-8), queda como sigue 

má + cá+ka=f (2.12-11) 


sistema de ecuaciones diferenciales formalmente análogo al (2.8-1); la única 
diferencia existente entre ambos estriba en el hecho de que las matrices m y c no 

tienen porqué ser diagonal. 

m y c son las denominadas matrices de masa y amortiguamiento 

consistente con la matriz de rigidez k , debido a que la discretización se consigue 

utilizando las mismas funciones de interpolación en ambos casos. 

Como ya se ha mencionado previamente la matriz de masa consistente o 
congruente es llena, lo que contrasta con las matrices de masa diagonales 
correspondientes a las masas concentradas. Debido a que es preferible trabajar con 
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este ultimo tipo de matrices, se suelen concentrar los valores obtenidos por algún 
procedimiento más o menos heurístico. Aunque ello introduce una inconsistencia en 
la formulación teórica tiene ventajas computacionales no despreciables y, dependiendo 
del tipo de problema, puede producir buenos resultados. 


Así, por ejemplo, a veces suele tomarse como matriz de masa la siguiente 


M = 


m l 0 ... 0 
0 m 2 ... 0 


0 


0 ... m 


N 


donde 


E m u 

<=1 


M 


M - masa total del sistema 


N 

E m a = traza de I a matriz consistente 

¿= i 


Una idea que se aplica con frecuencia y que es la más sencilla consiste en 
construir la matriz de masa del sistema suponiendo la masa concentrada en los g.d.l. 
de traslación; la masa asociada al giro se supondrá nula si considera la hipótesis de 
concentrar la masa en puntos sin inercia de rotación. Para ello se considera 
subdividida a la estructura en elementos; la masa de cada elemento se supone 
concentrada en los nudos extremos, asignando a cada nudo la mitad de la masa de 
elemento, tal y como se indica en el ejemplo 2.3. 


EJEMPLO 2.3 


Estudiar las vibraciones laterales de un edificio formado por una serie de 

pórticos espaciados 5 m. El peso de los forjados es de 500 y 250 kg/m 2 , y el de los 
muros 100 kg/m 2 . 

DATOS: * E = 2.1 x 10 6 kg/cm 2 

* Ij = 10200 cm 4 

* I 2 = 4400 cm 4 
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Solución : 


En la figura E.2.3-1 se presenta el esquema de uno de los pórticos interiores 
del edificio 



Fig. E.2.3-1 

* Peso de los forjados 

A p = 5 x 10 = 50 m 2 



= 12500 kg = 12.5 x 10 4 N 
m F2 = 12.5 x 10 3 kg 
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* Peso de los muros 



A m2 = 5 x 3 = 15 m 2 

P m2 = 100 x 15 tti 2 = 1500 kg = 15000 N 
m 2 

m m2 = 1500 k S 



A mI =5x5 = 25 m 2 

P ml = 100 -fe x 25 m 2 = 2500 kg = 25000 N 
m 2 

m ml = 2500 kg 


A la vista de lo anterior el modelo utilizado de masas concentradas sería el 
indicado en la figura E.2.3-2 



Fig. E.2.3-2 
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por lo que las masas concentradas en los grado de libertad de traslación son 


m i = ftl FI + m rn2 + tn m3 

= 25000 + 1500 + 2500 = 29000 kg 

m 2 = m n + m m2 = 

= 12500 + 1500 = 14000 kg 


• m 2 ■ m F2* m n,2 
t*1 ■ m F1* m m2* ra m1 


777/77 


Fig. E.2.3-3 


pudiendo escribirse la matriz de masa en la forma 


m 


m l 0 
0 m 2 



0 

14 


siendo la matriz de rigidez 






obtenida por el mismo procedimiento utilizado en el ejemplo E.2. 1. Así, para los dos 
grados de libertad indicados, la columna i son las fuerzas que aparecen en los g.d.l. 
j al dar un desplazamiento unidad en el g.d.l. i, manteniendo impedidos el resto de 
los grados de libertad. De este modo 
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- I a Columna : k¡, , j = 1 y 2 



- 2 a Columna : kjj , j = 1 y 2 

# 1 7^ * 1 tK k 22 



y como 

* L i = 5 m 

EI X = 2 10 6 x 10200 cm 4 = 204 x 10 8 kg cm 2 = 204 x 10 5 N m 2 
cmr 

* L 2 = 3 m 

EI 2 =2 10 6 x 4400 cm 4 = 88 * 10 8 kg cm 2 = 88 x 10 5 N m 2 
cm 2 
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resulta 


24 EL 94 y 204 v 10 5 N <■ N 

* k, = - = — - i-- — = 39.168 jc 10 5 — 

5 3 m m 


‘ I? 


= ^^2 = 24 x 88 £ 10 5 W = ?8 222 x 10 5 ü 
¿ 3* ™ ™ 


Las frecuencias propias se obtendrán al imponer la condición 

I k - w 2 m I = O 


que, en este caso, puede expresarse como 

(k t +k¿) - o 2 m l - k 2 


-K 


K ~ « 2 m 2 


= O 


lo que proporciona una ecuación de segundo grado en o? 2 

^k x +k^ - o) 2 ntjj jfcj - (o 2 tm 2 J - k% = O 


es decir 


2 [ m 2 (fcj + k¡) + m l k 2 ] ± \J[ m 2 ( k x + k 2 ) + m, ] 2 - 4 tn x m 2 

2 m l m 2 


para este caso particular se obtiene 

v>\ = 85.995 

w 2 = 877.529 

siendo entonces las frecuencias propias 

< 0 j = 9.273 rd/seg ; = 0.678 seg 

« 2 = 29.623 rdjseg ; T 2 = 0.212 seg 
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